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Preface 


Many  people  have  contributed  to  the  success  of  this  project  and  I  would  like  to  note  the  contributions 
of  the  individual  co-authors.  The  project  itself  is  a  testament  to  the  value  that  shorter  term  personnel 
can  add  to  a  program.  The  first  two  authors  were  not  full  time  employees  of  the  AFRL. 

•  Mark  Duell  came  to  us  in  the  summer  of  2005  as  a  “Space  Scholar”  (the  AFRLA^S  summer 
student  intern  program)  from  Embry-Riddle  Aeronautical  University,  between  his  junior  and 
senior  year.  His  job  the  first  year  was  to  assemble  and  test  the  ionic  anemometer  at  one 
atmosphere  and  in  a  bell  jar.  He  returned  after  his  senior  year  to  assemble  the  low  density  wind 
tunnel  we  used  to  test  and  to  perform  a  velocity  calibration  at  the  desired  altitude.  He  probably 
would  have  completed  that  if  the  laboratory  machine  shop  had  not  been  assigned  to  higher 
priority  tasking.  He  is  currently  employed  by  General  Electric  Infrastructure,  Aviation,  Eynn, 
MA. 

•  Larry  Saupe,  Maj  (S)  USAFR,  began  performing  his  reserve  duty  with  us  in  the  spring  of  2005. 
He  was  an  employee  of  General  Electric  Global  Research,  Niskayuna,  NY.  Larry  brought  an 
impressive  professional  background  in  aeronautical  and  mechanical  engineering  and  a  lot  of 
experience  with  problems  associated  with  power  generation,  heat  transfer,  and  fluid  flow.  In  the 
two  and  a  half  years  he  has  been  with  us,  serving  only  24  days  per  year,  Larry  has  contributed 
immensely  in  the  area  of  tunnel  configuration,  pump  and  fan  selection,  heat  transfer  analyses, 
tunnel  operation,  and  error  analysis. 

•  Brent  Barbeau,  Capt  USAF,  has  been  working  on  related  parts  of  this  project  since  his  arrival  in 
May,  2004.  He  originally  participated  in  the  selection  process  for  the  high  altitude  anemometer. 
Next  he  worked  with  Maj  Saupe,  and  mentored  Mark  Duell  through  his  two  summers  with  the 
project.  When  Mark  left,  Brent  worked  with  Larry  to  complete  the  tunnel,  and  he  ran  the  final 
testing  and  documentation  of  the  ionic  anemometer. 

•  Kris  Robinson  is  an  Electrical  Engineer  and  in-house  contractor  working  for  Utah  State 
University.  He  has  been  a  major  contributor  to  the  electrical  aspects  of  the  project,  to  name  some 
of  them:  initial  assembly  of  power  supply  components  to  the  anemometer,  testing  it  and  the  other 
candidate  instruments,  assembling  the  data  acquisition  system,  and  assisting  in  the  data  analysis. 
He  worked  together  with  Capt  Barbeau  running  the  final  checkout  experiments. 

There  were  also  many  other  people  who  contributed  to  the  program.  They  include,  Paul  Tracy, 
AFRLWSBY,  Lab  Manager  and  Section  Chief;  Will  Thom,  AFRLWSBX  who  lent  us  much  of  the 
equipment  needed  for  this  project,  and  who  gave  us  lots  of  great  advice;  Jack  Ignazio, 

AFRLWSBXR,  who  supplied  the  bell  jar  vacuum  assembly  and  assisted  us  in  its  use;  Radian  Belu, 
Electrical  Engineer  and  AFOSR  Summer  Research  Fellow,  currently  at  the  Desert  Research  Institute, 
NV,  who  assisted  with  many  high  voltage  aspects  of  the  problem;  Emmett  Bishop,  Electrical 
Engineering  Intern  from  Wentworth  Technical  Institute  for  circuit  design;  and  John  Redwood, 
machinist  with  the  AFRLWS  machine  shop,  who  made  the  ionic  anemometer.  Lts.  Mike  Walker  and 
James  Baird  were  active  participants  in  the  earliest  phases  of  the  program.  We  are  also  indebted  to 
Dr.  Arje  Nachman,  our  AFOSR  sponsor  for  this  basic  research  project. 


1.  INTRODUCTION 


Aircraft  occasionally  encounter  turbulence  in  flight.  There  has  been  considerable 
effort  to  try  to  foreeast  the  occurrence  of  this  turbulenee  by  NASA,  the  FAA,  NOAA,  and 
NCAR.  Systems  are  in  plaee  and  are  being  evaluated  to  determine  the  best  ways  to 
detect  turbulence  and  warn  pilots'.  High  altitude  aireraft,  sueh  as  the  U-2  also  encounter 
turbulenee,  whieh  is  often  Clear  Air  Turbulence  or  CAT,  not  assoeiated  with  any 
thunderstorms  or  other  convective  activity.  Often  the  turbulenee  appears  to  be  assoeiated 
with  flight  over  mountains.  It  appears  that  a  major  cause  is  the  gravity  or  buoyancy 
waves  generated  by  flow  over  the  mountains,  whieh  grow  in  amplitude  as  they  ascend, 
until  they  are  either  absorbed  into  the  flow  by  a  critieal  layer,  or  until  they  break  into 
turbulence^.  The  exact  nature  of  the  wave  breaking,  and  the  identification  of  exactly 
what  conditions  are  dangerous  to  aireraft  are  not  fully  understood. 

Most  turbulence  foreeast  models  are  not  set  up  to  operate  at  these  altitudes,  well  above 
commercial  airline  flight.  It  is  not  eertain  that  operating  the  models  at  the  higher  altitude 
would  in  fact  capture  these  events,  due  to  resolution  at  the  high  altitudes  and  other 
numerical  problems.  There  are  only  a  small  number  of  foreeast  tools  available  to  the 
high  altitude  pilot.  One  is  the  Mountain  Wave  Foreeast  Model  (MWFM)  operated  by  the 
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Naval  Researeh  Laboratory  ,  which  uses  ray  traeing  methods  to  project  mountain  waves 
through  global  foreeast  models.  Other  eandidate  foreeast  models  will  be  examined  to 
determine  their  ability  to  detect  CAT  conditions  at  high  altitudes.  There  is  always  a 
problem  of  validating  the  output  of  any  turbulence  model  (see  Ref  1).  Ideally  the  output 
of  the  model  should  be  compared  to  aetual  data,  but  when  turbulenee  is  forecast,  no  one 
wants  to  fly  through  it. 

Recently"',  it  was  proposed  that  the  US  should  investigate  the  utility  of  operating  a 
blimp-like  object  at  altitudes  above  60  000  ft.  This  idea  developed  into  the  High  Altitude 
Airship  (HAA),  an  Advanced  Concept  Teehnology  Demonstration  project  managed  by 
the  U.S.  Army,  funded  through  the  Missile  Defense  Agency.  Originally  slated  to  fly  in 
2006,  it  has  encountered  some  development  problems,  and  is  currently  scheduled  to  fly  in 
2008.  The  HAA  coneept  is  to  ferry  itself  from  a  home  base  to  an  operating  loeation,  then 
maintain  station  keeping  for  long  periods  of  time.  In  studies  of  wind  elimatology  at  the 
target,  it  has  been  estimated  that  the  airship  should  be  able  to  sustained  headwinds  up  to 
30  knots.  While  considering  the  problems  associated  with  the  validation  of  high  altitude 
CAT  foreeasts  and  the  general  lack  of  long  term  observations  at  high  altitude,  the  HAA 
appears  to  be  an  ideal  platform  to  study  waves  and  turbulenee  at  high  altitudes.  Not  only 
would  this  help  in  the  verification  of  wave  and  turbulenee  forecast  models,  but  it  would 
also  provide  the  environmental  specifieation  for  future  stratospheric  airship  payloads 

2.  REQUIREMENTS 

The  design  eondition  will  be  for  a  standard  atmosphere  at  60  000  ft  (18  288  m): 
T=390  R  (217  K),  P=1.05  psi  (7,230  Pa),  p  =0.000226  slugs/ft^  (0.116  kgW).  The 
mean  speed  of  the  HAA  based  on  elimatology  requires  a  sustained  cruise  speed  of  30 
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knots^  or  50  ft/s  (15  m/s).  The  desired  minimum  wind  varianee  is  3ft/s  (Im/s).  Minimal 
temporal  sensitivity  would  be  to  deteet  a  flow  disturbanee  down  to  an  eddy  size  of  1  m  at 
30  m/s,  or  0.033  s.  This  would  require  at  least  60  Hz  for  two  data  points  per  eddy  at  the 
maximum  cruise  speed.  The  reference  altitude  is  60  000  ft. 

3.  ALTERNATIVES 

Four  alternatives  were  considered  for  use  as  a  high  altitude  velocity  turbulence  sensor: 
the  pitot-static  tube,  the  hot  wire  anemometer,  the  propeller  anemometer,  and  the  ionic 
anemometer,  which  was  used  by  French  scientist,  Jean  Barat  for  many  high  altitude 
measurements.  The  suitability  of  each  for  the  required  mission  is  analyzed  below. 

3.1  Pitot  -  Static  Tube 

A  pitot-static  tube  is  a  flow  velocity  meter  which  is  capable  of  measuring  fluid 
velocities  as  a  localized  point  (as  opposed  to  an  averaged  velocity  across  a  larger 
section).  The  pitot-static  tube  senses  pressure  at  two  locations:  one  in  front  -  ideally  the 
stagnation  point  where  the  flow  stops  yielding  the  “stagnation”  or  “total”  pressure,  and 
one  on  the  side  -  ideally  at  a  spot  where  the  pressure  is  indicative  of  “free  stream”  or 
“static”  pressure.  For  lower  speed  flows,  the  important  quantity  is  the  difference  in  the 
two  pressures,  known  as  the  “dynamic  pressure”,  which  is  typically  measured  with  a 
differential  pressure  transducer.  The  fluid  velocity  can  be  obtained  from  the  difference  in 
pressure  or  dynamic  pressure  using  incompressible  fluid  theory.  These  probes  are  also 
used  for  high  speed  flows,  where  Mach  Number  can  be  deduced  from  the  ratio  of  the 
measured  pressures,  and  speed  is  determined  from  a  speed  of  sound  calculation  that 
requires  the  air  temperature. 

A  flow  can  be  considered  incompressible  if  its  velocity  is  less  than  30%  of  its  sonic 
velocity.  For  such  a  fluid,  the  Bernoulli  equation  describes  the  relationship  between  the 
velocity  and  pressure  along  a  streamline: 


^  +  z  +  ^  =  C  (1) 

ig  pg 

Evaluated  at  two  different  points  along  a  streamline,  the  Bernoulli  equation  yields: 

^  +  z  +^  =  ^  +  z  m 

2g  ^  pg  2g^  2  ^  pg  (2) 


If  zi  =  Z2  and  point  2  is  a  stagnation  point,  i.e.,  V2  =  0,  the  above  equation  reduces  to: 


2^  _l_  _Pi_  _ 

1  p  p 


(3) 
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For  instrument  sensitivity,  we  eompute  the  pressure  difference  for  a  3  ft/s  (~1  m/s) 
variance  with  various  mean  wind  speeds  up  to  50  ft/s  (~15  m/s).  Solving  for  the  pressure 
difference; 

^P  =  ^p[v+v')^  (4) 

The  results  are  plotted  in  Figure  1 .  For  the  maximum  mean  wind  velocity,  the 
difference  in  pressure  is  a  bit  over  0.002  psi  (17Pa),  which  decreases  as  wind  velocity 
decreases.  While  pressure  differences  of  this  low  magnitude  are  measurable,  they  are 
below  what  is  commonly  available  for  operational  -  long  endurance  use. 

3.2  Hot  Wire  Anemometer 

Since  the  main  purpose  of  the  project  is  to  measure  turbulence,  a  common  laboratory 
solution  is  the  hot  wire  anemometer,  which  was  briefly  considered  for  this  application.  A 
hot  wire  does  lose  some  performance  as  air  density  decreases  and  the  heat  transfer 
coefficient  is  reduced^.  Use  at  FIAA  altitudes  reduces  the  frequency  response  to  about 
one  half  of  the  sea  level  response.  The  primary  concern  is  the  need  for  long  endurance. 
Hot  wires  must  be  very  thin,  typically  1pm,  to  achieve  the  required  frequency  response 
for  turbulence.  The  net  result  is  a  rather  fragile  instrument,  that  would  probably  not  be 
suitable  for  long  endurance.  The  comer  frequency  (the  highest  frequency  of  good  data) 
of  a  hot  wire  for  a  range  of  atmospheric  density  ratios,  up  to  1  (sea  level  density),  and  for 
a  range  of  wire  diameters  from  5  to  50  pm  is  shown  in  Figure  2.  The  approximate 
density  at  HAA  altitude  is  the  at  the  left  side  of  the  curve. 

The  primary  cause  of  decreased  frequency  response  is  the  increase  in  wire  diameter. 
In  order  to  meet  the  60  Hz  requirements,  the  wire  must  be  no  greater  than  7.5  pm  in 
diameter,  which  is  very  thin  from  an  endurance  standpoint. 


Mean  Velocity  (ft/s) 


Figure  1.  Pressure  difference  for  a  pitot- 
static  probe  at  HAA  altitudes  for  wind 
velocities  up  to  5  Oft/s  for  a  velocity 
difference  of  3ft/s. 
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Figure  2.  Corner  frequency  of  a  hot  wire 
anemometer  vs.  atmospheric  density  ratio  for 
various  wire  diameters 
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3.3  Propeller  Anemometer 

Anemometer  devices  have  existed  in  various  forms  over  the  past  several  centuries. 
One  of  the  more  popular  recent  designs  is  the  propeller  anemometer,  a  device  consisting 
of  several  propeller  blades  that  rotate  in  the  presence  of  an  airflow.  An  example  of  one 
such  device  is  depicted  in  Figure  3.  The  rate  at  which  the  propeller  spins  when 
encountering  an  airflow,  is  used  to  quantify  the  magnitude  of  that  airflow. 


Figure  3.  Cross  sectional  view  of  a  propeller  through  one  of  the  blades  at  a  located  a 
distance  r  from  the  axis  of  rotation.  Lift,  L,  and  drag,  D,  are  referenced  to  the  zero  lift  line. 
If  the  propeller  is  spinning  at  just  the  right  speed  to  cause  no  thrust  or  drag  on  the  propeller 
shaft,  then  all  the  airflow  will  be  along  the  axis  of  the  blade.  Vp  is  the  air  velocity 
approaching  the  propeller,  V  is  the  rotational  velocity  approaching  the  blade,  Vt  is  the  total 
velocity  flowing  over  the  blade  element. 

3.3.1.  Propeller  Anemometer  Theory 

At  moderate  altitudes,  a  lightweight  propeller  anemometer  with  a  high  number  of 
blades  exhibiting  a  large  surface  area  can  measure  airflows  of  small  intensity,  or  small 
fluctuations  in  the  mean  airflow.  The  speed  of  the  airflow  is  obtained  by  measuring  the 
rate  at  which  the  airflow  causes  the  propeller  to  rotate.  A  propeller  anemometer  is  a 
logical  candidate  for  the  HAA;  however  the  friction  and  inertia  of  the  overall  system,  and 
the  low  air  density  at  the  desired  altitude  and  low  airspeed  could  make  it  an  ineffective 
device  on  the  HAA. 
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To  predict  the  performance  of  a  typical  propeller  anemometer  on  the  HAA,  a  simulation 
of  the  propeller  was  examined  using  MATLAB®^.  The  model  was  run  with  the 
parameters  of  a  commercial  high  quality  turbulence  sensing  propeller  anemometer  made 
by  the  R.M.  Young  Company.  Results  of  the  simulations  were  compared  to  wind  tunnel 
test  results  to  fine-tune  the  model. 

3.3.2.  Propeller  Anemometer  Model 

To  derive  a  mathematical  model  for  a  propeller  anemometer,  we  begin  with  the 
conservation  of  angular  momentum  equation,  the  rotational  analog  of  Newton’s  2"^^  Law: 

(5) 

where  J  is  the  total  moment  of  inertia  of  the  system,  Q’  is  the  angular  acceleration  and 
EM  is  the  sum  of  the  moments  on  the  propeller  system,  all  about  the  axis  of  rotation. 

An  anemometer  has  three  contributions  to  the  inertia:  propeller,  shaft,  and  the  motor. 
The  moments  on  the  anemometer  are  the  torque  from  the  propeller,  T,  and  any  rotational 
friction,  expressed  as  a  constant,  bj,  multiplied  by  the  rate  of  rotation: 

{Jp+Js+Jm)^'=T -b^Q.  (6) 

Since  our  main  purpose  is  to  see  if  an  anemometer  is  even  feasible  at  the  altitude  at 
which  the  HAA  will  fly,  we  first  make  some  simplifying  assumptions.  If  an  ideal, 
frictionless  anemometer  does  not  give  us  the  needed  results,  then  it  will  not  even  be 
necessary  to  try  computing  one  with  all  of  the  variables  added.  Therefore,  we  simplify 
the  equation  by  ignoring  damping  and  the  inertias  of  the  shaft  and  the  motor,  assumed  to 
be  much  less  than  the  propeller. 


(7) 

o 

Following  blade  element  theory  ,  for  a  thin  section  of  the  propeller  blade  at  a  distance 
r  from  the  axis  of  rotation,  the  equation  for  incremental  torque  on  a  propeller  is: 

dT  =  r\dLsm{(l))  +  dDcos{(l))\  (8) 

where  dL  is  the  incremental  lift,  andfiD  is  the  incremental  drag,  which  will  be  ignored 
for  this  study.  Lift  is  perpendicular  to  the  zero  lift  line  through  the  blade  element,  and  (j) 
is  the  angle  between  the  lift  vector  and  the  thrust  vector,  which  is  also  the  difference 
between  the  local  blade  pitch  angle  {9)  and  the  angle  of  attack  (a).  Using  the  standard 
equation  for  lift  in  terms  of  a  coefficient  of  lift,  where  is  the  local  flow  velocity,  is 

the  lift  coefficient,  c  is  the  chord,  b  is  the  number  of  blades,  and  dr  is  the  increment  of 
radius: 


dT  =  ^pV^bc[C^^sm{(j))\rdr 


(9) 
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Since  is  assumed  to  be  a  linear  function  of  a  ,  it  can  be  substituted  using  the 
following  equation: 


C^=a{e-(p')  (10) 

In  this  equation,  a  represents  the  lift  curve  slope,  assumed  eonstant.  The  equation  now 
becomes: 


dT  =  pV^bc[a{6  -  (p)  sin(^zi)]rt/r 


(11) 


As  mentioned  earlier,  is  the  section  flow  velocity,  which  is  the  vector  sum  of  the 
axial  flow  into  the  propeller  (Vp),  and  the  angular  flow  veloeity  (V). 

V,  =  V^+V  (12) 


The  angular  flow  veloeity  is  the  produet  of  the  angular  veloeity  of  the  blade  and  the 
radius  at  that  point  of  the  blade,  V  =  Qr  .  The  angle  between  Vp  and  V  is  .  Using 
trigonometry  we  arrive  at  the  following: 


tan(^zi)  =  ^ 
Llr 


V 

sin(^zi)  =  — 


Assuming  is  a  small  angle  gives  the  following  results:  tan(^zi)  =  <!> ,  (j)  =  V  ^  /Q. 


and  «  F  «  Qr .  Note  that  this  assumption  can  not  be  used  when  the  blades  have 

stopped  or  are  rotating  slowly  relative  to  the  wind  speed.  The  Torque  equation  now 
becomes: 


1 


dT  =  —  c£l  r  abc 
2 


{ 


Qr 


A 


F„ 


Qr 


rdr 


(13) 


The  equation  for  the  pitch  of  an  ideally  twisted  blade  is:  6  =  — — ,  where  6^  is  the  pitch 

r 

angle  at  the  blade  tip,  and  R  is  the  entire  radius  of  the  blade.  Substituting  for  6  and 
multiplying  the  terms  yields: 


dT  =  pD.rabcV -^prabcVp)dr  (14) 

The  full  propeller  torque  is  the  integral  of  the  incremental  torque  from  the  inner  radius, 
assumed  to  be  zero  to  R,  the  full  blade  radius: 

R  R 

T  =  \ pQ.abcVpO,R^rdr  -  \ pabcVj' ^rdr  (15) 

0  0 


After  the  integration,  and  faetoring  the  equation  becomes: 
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(16) 


1  ,  Vp 

T  ^-pabcVB,R\Q. - 

The  differential  equation  for  rotational  veloeity,  Q  ,  is  now: 

1  V 

JQ' — pabcVe,R\D. - ^)  =  0  (17) 

This  is  a  nonlinear  differential  equation.  The  term  V^^jO^R  is  the  rotational  veloeity  that 

results  in  no  thrust.  The  differenee  between  aetual  rotation  and  the  no  thrust  rotation  is 
the  “error”  term  that  eauses  the  blade  rotational  speed  to  inerease  or  deerease.  Also  note 
that  Vp  and  p  are  terms  in  the  gain  that  are  available  to  eorreet  the  speed  error  by 
ehanging  the  rotation  rate  of  the  propeller.  MATLAB®  was  used  to  solve  the  equation. 

3.3.3.  Propeller  Anemometer  Tests: 

In  order  to  ealibrate  the  model,  the  simulation  was  eompared  to  aetual  results  from  a 
wind  tunnel  experiment.  The  anemometer  used  for  the  study  was  one  leg  of  a  triple 
propeller  “Gill  UVW  Anemometer”  sold  by  R.M.  Young  Company^,  Model  27005.  The 
propeller  is  a  four  blade  expanded  polystyrene  propeller,  model  08274  EPS  Propeller, 
made  by  the  R.M.  Young  Company.  The  properties  of  this  propeller  were  input  as  the 
variables  in  the  MATLAB®  model,  and  a  simulation  was  run  using  the  same  initial 
eonditions  as  the  wind  tunnel  testing.  These  variables  are  shown  below.  The  parameters 
of  the  blade  are:  R  =  0. 1  Im,  e  =  0.05m,  mass  =  0.01 1kg.  The  propeller  inertia  depends 
on  the  mass  and  radius  of  gyration  of  the  propeller  (Jp=mRg^).  For  the  radius  of  gyration, 
we  used  45%  of  the  radius.  The  piteh  angle  was  listed  as  0.425rad  in  the  eatalog.  The 
slope  of  the  lift  eurves,  a,  is  assumed  to  be  6.8/rad,  a  representative  value  for  propellers. 
The  test  was  run  at  approximately  sea  level  eonditions,  with  p  =  \  .2kg/m  ,  with  a  wind 
speed  of  7.6m/s.  For  these  eonditions,  the  piteh  angle  had  to  be  adjusted  to  0.51  radians 
to  aehieve  the  eorreet  final  speed. 

A  eomparison  of  the  wind  tunnel  data  to  the  model  is  shown  in  Figure  4.  The  results, 
using  the  above  variables  shown  as  the  blue  dashes,  labeled  “Model  1”  are  somewhat 
disappointing.  In  talking  to  the  manufaeture,  we  found  that  some  of  our  original 
assumptions  were  quite  optimistie.  For  one  thing,  the  eombined  moment  of  inertia  is 
apparently  more  than  what  was  used.  Also,  the  eompany  eites  a  non-negligible  frietion 
when  the  propeller  is  stopped.  The  eompany  eites  a  sensitivity,  or  distanee  eonstant.  as 
1.0  m;  whieh  is  an  indieation  of  the  time  to  get  to  63%  of  full  veloeity  eompared  to  the 
time  it  takes  for  1  m  of  fluid  to  pass  the  blade.  At  this  veloeity,  the  time  should  have 
been  0.13  s.  The  moment  of  inertia  for  the  blade  was  inereased  to  aehieve  this  result, 
whieh  is  also  shown  on  the  plot  as  Model  2.  This  is  still  about  twiee  as  fast  as  the  aetual 
data,  but  sinee  we  should  have  started  at  something  other  than  zero  speed,  and  sinee  we 
had  an  older  anemometer,  we  deeided  to  use  the  eompany  provided  sensitivity. 

Next  the  simulation  was  run  for  the  design  altitude,  by  ehanging  the  density.  The 
simulation  was  run  for  a  Im/s  step  inerease  in  wind-speed  starting  from  0,  10,  20,  and 
30  m/s.  The  times  to  the  63%  point  for  sensitivity  ealeulations  were  10  s,  1  s,  0.45  s,  and 
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Figure  4.  Propeller  rotation  rate  vs  time. 
Comparison  of  simulations  to  wind  tunnel  data. 
Wind  tunnel  data  is  represented  by  the  dots, 
Original  model  is  dashed.  The  solid  curve  uses 
manufacturers  sensitivity. 
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Figure  5.  Estimated  time  to  reach  63%  of 
steady  state  vs.  wind  speed  (velocity)  at  65  000 
ft  altitude,  using  the  manufacturer’s  sensitivity 
in  the  simulation. 


0.4  s  respectively,  shown  in  the  plot  in  Figure  5.  This  equates  to  a  distance  constant  of 
around  10  m,  an  order  of  magnitude  higher  than  the  sea  level  value  quoted  by  the 
manufacturer.  Clearly  the  low  density  takes  a  toll  on  the  responsiveness  of  the  system, 
and  would  prohibit  the  instrument  from  meeting  the  desired  response  rate  of  the  velocity 
turbulence  sensor.  The  use  of  propellers,  shafts,  and  sensors  with  lower  moments  of 
inertia  is  possible,  but  this  action  would  impact  instrument  durability  and  serviceability. 


3.4.  Ionic  Anemometer 


A  one  dimensional  ionic  anemometer,  that  is,  capable  of  measuring  only  one  velocity 
component,  consists  of  three  parallel  wires;  one  emitter  wire  and  two  collector  wires  on 
opposite  sides  of  the  emitter  at  equal  distance.  When  a  large  potential  difference  (voltage) 
is  applied  between  the  emitter  and  collector,  some  air  particles  near  the  emitter  are 
ionized  (believed  to  be  primarily  (Fl20)n02^  and  (Fl20)nN2^).  This  voltage  is  the 
activation  voltage;  it  is  on  the  order  of  tens  of  kilovolts  for  a  1cm  separation  at  1 
atmosphere.  The  electric  field  between  the  emitter  and  collectors  accelerates  the  ions 
toward  the  collectors,  producing  a  small  current  flow  (on  the  order  of  microamperes). 

The  acceleration  of  these  ions  occurs  in  negligible  distance  and  they  travel  the  majority  of 
the  distance  at  a  steady  state  velocity,  Vi,  determined  by  the  strength  of  the  electric  field, 
density  of  the  air,  and  humidity  of  the  air.  One  must  insure  that  the  voltage  stays  low 
enough  to  avoid  arcing. 

3.4.1  Ionic  Anemometer  Theory 

In  reported  high  altitude  experiments,  J.  Barat^°  used  two  ionic  anemometer  designs  to 
analyze  the  stratospheric  winds;  one  used  pulsed  voltage  to  create  ions,  the  second  used 
continuous  ionic  flow.  We  focused  on  the  latter.  He  found  that  the  design  operated  well 
between  200  and  10  hPa  (40  000  to  100  000  ft),  had  a  detection  limit  better  than  1  cm  s'\ 
an  overall  accuracy  better  than  10%,  a  response  time  constant  lower  than  10'  s,  and  was 
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sufficient  to  measure  CAT  in  the  stratosphere.  Barat  gives  the  equation  for  ion  veloeity, 

Vi. 

(18) 

where  is  the  mobility  of  the  ions  and  E  is  the  strength  of  the  eleetrie  field.  The 
mobility  fj.^  is  a  funetion  of  density: 


where  /u^  is  the  redueed  mobility,  a  funetion  of  humidity,  but  independent  of  dry  air 
density,  p  .  As  the  air  density  deereases,  the  mobility  of  the  ions  inereases,  and  thus  the 
required  E  field  strength  for  a  desired  veloeity  deereases. 

Consider  the  one  dimensional  anemometer  in  Figure  6  (designed  to  measure  veloeity 
in  only  one  direetion).  Ions  are  released  from  the  eentral  emitter  in  all  direetions,  with  a 
slow  initial  veloeity  that  rapidly  eomes  up  to  the  ion  veloeity  relative  to  the  air  veloeity. 

If  the  same  number  of  positive  ions  went  to  the  left  as  to  the  right,  the  ion  density  would 
be  higher  on  the  upstream  side,  whieh  would  loeally  modify  the  eleetrie  field.  This  must 
eneourage  more  ions  to  drift  downstream,  tending  to  equalize  the  ion  density.  The 
density  of  the  positive  ions  must  be  approximately  eonstant  within  the  instrument,  sinee  it 
has  been  found  that  for  the  eolleetors  upstream  (1)  and  downstream  (2): 

h  ^V.-V^  and  (20) 

Assuming  the  eonstants  are  the  same  (symmetrieal  instrument),  this  yields: 

Ei^(V,-VJ-(V,+VJ  =  2V^ 

(21) 

^i^(Vi-VJ  +  (V,+VJ  =  2V^. 

(22) 

Dividing  Eq.  (21)  by  Eq.  (22)  and  rearranging 
yields: 

(23) 

where  K  is  the  sensitivity  eonstant.  In  normal  use, 

V.  should  be  eonstant  in  the  presenee  of  ehanging 
wind  eonditions. 
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Figure  6.  Sehematie  of  a  one 
dimensional  ionie  anemometer. 


As  a  preliminary  test  we  need  to  find  the  activation  voltage  and  size  of  the  region  that 
is  usable  for  measurements.  For  sensitivity  calibration,  we  need  to  find  the  slopes 
(sensitivity  constant)  and  intercepts  (due  to  non-symmetry)  of  velocity  vs.  AZ/I/  curves 
for  several  voltages.  For  density  calibration,  we  needed  to  find  the  activation  voltage  and 
verify  operation  in  low  density  humid  and  dry  air.  The  transit  time  for  each  of  the  ions 


Figure  7.  Top  and  cross-sectional  view  of  the  ionic 
anemometer.  Units  are  in  inches.  Center  electrode  is 
tensioned  by  spring  in  the  bottom  housing  and  friction  from 
the  machine  screw  on  the  top  housing. 


Figure  8.  Photograph  of  the 
ionic  anemometer  attached  to 
the  high  voltage  leads  and 
ammeters. 


from  emitter  to  collector  is  believed  to  be  very  short,  giving  the  instrument  a  rapid 
response  time  to  changing  wind  conditions  allowing  for  the  desired  frequency  of  response 
(400Hz).  The  von  Karman  frequency  of  the  flow  passing  over  the  collectors  ranges  from 
6  Hz  at  1  ft/sec  to  2  kHz  at  70  ft/sec  and  may  cause  fluctuations  at  some  velocities  in  the 
desired  usable  range.  Adding  additional  collectors  (e.g.  to  determine  velocity  in  the 
perpendicular  plane)  increases  the  strength  of  the  electric  field,  which  changes  the  ionic 
drift  speed,  activation  voltage,  and  current. 


3.4.2.  Ionic  Anemometer  Hardware 


The  overall  shape  and  design  of  the  instrument,  shown  in  Figures  7  and  8,  were  based 
on  Barat’s  design.  The  overall  length  of  the  instrument  is  0.15  m,  with  a  diameter  of  0.03 
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m.  It  is  desirable  to  have  a  very  fine  emitter  wire  (20  -  120  [am  is  typieal);  100  pm 
Tungsten  was  ehosen  for  durability  and  availability.  The  eolleetor  rods  are  struetural 
members,  and  thus  must  be  sized  to  withstand  antieipated  loads;  however  large  eolleetors 
ereate  higher  magnitude,  lower  frequeney  turbulenee,  so  minimizing  their  size  is  desired. 
Aluminum  rods  1.6  mm  in  diameter  were  seleeted  for  the  eolleetors  sinee  they  are 
expeeted  to  be  suffieient  and  available.  The  overall  length  of  the  eolleetors  is  0.14  m, 
with  an  exposed  length  of  0.12  m.  The  emitter  wire  has  an  exposed  length  of  0.09  m; 
there  is  some  shielding  at  either  end  so  that  slightly  off-axis  winds  do  not  affeet  total 
eurrent  readings.  Both  the  emitter  and  eolleetors  are  free  from  any  surfaee  eoatings  or 
finishes  in  order  to  maximize  eleetrie  eonduetivity.  Some  ionie  anemometers  use 
radioaetive  eoatings  on  the  emitter  to  aid  in  the  produetion  of  ions.  These  eoatings  greatly 
reduee  the  voltage  requirement  (by  approximately  two  orders  of  magnitude),  but  they 
also  limit  instrument  life  and  add  signifieant  eomplexity  to  instrument  eonstruetion, 
availability,  and  handling.  We  eleeted  to  use  higher  voltages  instead  of  radioaetive 
eoatings,  whieh  present  their  own  eoneerns.  When  dealing  with  high  voltage  (espeeially 
during  the  1  atmosphere  tests  on  the  ground),  suitable  spaeing  of  all  parts,  insulation,  and 
eonneetors  as  well  as  proper  safety  proeedures  must  be  used.  The  struetural  disks  at  eaeh 
end  must  be  eleetrieally  insulating  between  the  emitter  and  eolleetors.  Teflon  was 
seleeted  for  its  durability,  low  temperature  behavior,  availability,  and  laek  of  outgassing 
at  low  pressures.  Outgassing  was  also  a  eoneern  in  the  ehoiee  of  adhesives,  and  the 
outgassing  from  the  superglue  we  seleeted  was  verified  to  be  suffieiently  low  at  expeeted 
flight  pressures.  The  emitter  wire  is  tensioned  by  means  of  a  plate  soldered  to  the  end  of 
it  and  a  eompressed  spring  in  the  base  of  the  instrument.  This  forees  the  emitter  wire  to 
remain  straight  even  when  temperature  ehanges  eause  the  length  of  the  eolleetors  to 
ehange.  At  the  top  of  the  instrument  the  wire  is  held  in  plaee  by  frietion  between  a  serew 
and  the  Teflon  disk.  An  eleetrie  eonneetion  is  maintained  by  a  thin  eopper  shim  that  fits 
around  the  spring  and  provides  a  lead  to  the  outside. 

3.4.3.  Ionic  Anemometer  Preliminary  Testing 

The  anemometer  was  plaeed  in  a  test  stand  with  positive  high  voltage  eonneeted  to  the 
emitter  and  two  eolleetors  on  opposite  sides  eonneeted  to  separate  miero  ammeters.  The 
two  other  eolleetors  were  eonneeted  to  ground.  To  determine  the  aetivation  voltage,  the 
voltage  was  inereased  from  zero  until  one  of  the  ammeters  read  0.2  pa.  This 
voltage/eurrent  data  were  reeorded.  To  determine  the  eharge  eharaeteristie,  shown  in 
Figure  9,  the  voltage  was  further  inereased  to  125%  of  the  aetivation  voltage  and  the 
voltage/eurrent  data  were  reeorded  at  regular  intervals.  The  voltage  was  then  deereased 
until  shut-off,  reeording  the  data  at  regular  intervals. 

For  veloeity  ealibration  at  sea  level,  the  voltage  was  set  to  105%  of  aetivation,  and  the 
eurrents  were  reeorded  at  various  wind  speeds  (ineluding  no  wind).  This  was  repeated  at 
110%,  115%,  and  120%  of  aetivation  voltage.  The  eurrent  differenee  divided  by  the  sum 
was  plotted  versus  veloeity  in  Figure  10.  At  either  of  the  two  applied  voltages,  the 
response  to  veloeity  appears  quite  linear,  indieating  a  eonstant  sensitivity. 

Next,  the  probe  anemometer  was  tested  at  lower  pressures,  near  those  expeeted  at 
operating  altitude  near  60  000  ft.  The  instrument  was  plaeed  in  a  bell  jar,  and  pressure 
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was  drawn  down  to  73  hPa,  approximately  7%  of  standard  sea-level  atmospheric 
pressure.  First  the  single  collector  current  was  measured  as  a  function  of  voltage,  now 
much  lower,  as  expected.  The  pressure  in  the  jar  was  then  drawn  to  near  vacuum,  then 
refilled  to  73  hPa  with  dry  nitrogen  and  the  test  was  repeated.  The  results,  shown  in 
Figure  11,  were  similar  to  the  sea-level  curves,  but  at  significantly  lower  voltages.  While 
a  calibrated  flow  source  was  not  available  in  the  bell  jar,  a  computer  cooling  fan  that 
could  be  run  at  two  input  voltages  was  placed  near  the  anemometer.  The  current 
difference  over  the  sum  is  shown  for  4  applied  voltages  in  Figure  12. 


Figure  9.  Current  (pA)  vs.  applied  voltage  (kV)  for  a  single  collector  at  1010  hPa. 
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Figure  10.  Velocity  calibration  for  two  collectors  at  1010  hPa  and  two  applied  voltages, 

5.3  and  6.0kV. 

To  determine  the  effect  of  density,  the  activation  voltage  was  measured  as  described 
above,  at  densities  from  pressures  at  60  000  ft  down  to  those  at  around  20  000  ft.  Those 
results  are  shown  in  Figure  13.  The  activation  voltage  increased  from  IkV  to  3kV  as 
pressure  increased,  as  expected. 
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♦  Humid  Air 
■  Dry  Nitrogen 

I _ - 1  I 

Figure  1 1 .  Current  (fxA)  vs.  applied  voltage  (kV)  for  two  collectors  in  dry  N2  and 
humid  air  at  73  hPa. 


Figure  12.  Velocity  calibration  for  two  collectors  in  dry  N2  at  73  hPa  as  shown  at 
various  applied  voltages  to  the  emitter. 
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4.  LOW  DENSITY  TUNNEL  SELECTION  PROCESS 


4.1  Introduction 

As  part  of  the  High  Altitude  Airship  (HAA)  payload  program,  a  test  eell  is  desired  to 
replieate  several  in-flight  atmospherie  eonditions  that  may  be  experieneed  by  proposed 
test  instrumentation  during  flight.  This  test  eell  would  primarily  be  used  to  validate 
predieted  instrument  response  to  high  altitude,  low  density  eonditions. 

Key  parameters  for  this  test  eell  are  duplieation  of  the  air  density  and  veloeity 
eneountered  in  aetual  flight  eonditions.  Air  density  is  ehosen  to  be  that  whieh  would  be 
experieneed  during  Standard  Atmosphere^^  eonditions  at  60  000  ft.  Chosen  flight 
veloeity  is  30  knots  (whieh  has  been  assumed  as  30  knots  true  airspeed).  The  test  eell  has 
been  initially  desired  to  have  a  eireular  eross  seetion  with  an  internal  diameter  of  8.0 
inehes. 

Three  different  approaehes  are  outlined  in  this  report;  Once  Through,  Closed  Loop 
Vacuum  Pump  and  Closed  Loop  Fan.  Of  the  three.  Closed  Loop  Fan  appears  to  be 
initially  most  attraetive  due  to  its  simple  eonstruetion,  ease  of  use,  low  initial  eost  and 
low  eost  of  operation.  However,  depending  upon  equipment  availability  and  duration  of 
testing,  other  options  may  ultimately  prove  more  desirable.  The  three  eonfigurations  are 
shown  in  Figures  17,  18,  and  19. 


4.2  Test  Cell  Parameters 

Table  1  outlines  the  aetual  and  desired  test  eonditions  for  this  test  eell  (utilizing  air  as  a 
test  fluid). 

Table  1.  Aetual  and  desired  test  eell  eonditions 


Parameter 

StdAtm  @  60,000  ft 

Test  Cell  w/Air 

Pressure  (millibar)  mmHgA 

(72.312)  54.24 

(98.22)  73.67 

Temperature  (°C)  °F 

(-56.5)  -69.7 

(21.1)70.0 

Density  (kg/m^)  Ib/ft^  * 

(0.11628)  0.007262 

(0.11628)  0.007262 

Veloeity  (knots)  ft/s 

(30.0)  50.6 

(30.0)  50.6 

C  (knots)  ft/s 

(573.8)  968.5 

(668.7)  1128.7 

Maeh  # 

0.05228 

0.04863 

Vdot  8”  dia  duet  (m^/min)  ft^/min 

(30.03)  1060.5 

(30.03)  1060.5 

Mdot  8”  dia  duet  (kg/hr)  Ib/hr* 

(209.6)  462.1 

(209.6)  462.1 

Dynamie  AP  head  (cm-H20)  in-H20 

(0.1414)  0.05567 

(0.1414)  0.05567 

Reynolds  Number  (duet  diameter) 

25,652 

20,051 

*  Assumes  G  as  32. 1 74 ft/s^ 
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4.3  Once  Through  configuration 

A  Once  Through  configuration  is  initially  proposed  due  to  its  inherent  simple 
construction  and  a  readily  available  high  vacuum  source  (Busch  vacuum  pumps)  at 
Hanscom  AFRL. 

A  proposed  flow  schematic  is  provided  in  Figure  17.  For  this  configuration, 
atmospheric  air  is  drawn  into  the  system  through  one  or  more  critical  flow  nozzles. 
These  flow  nozzles  will  provide  two  important  functions;  1)  dramatically  lower  the 
stagnation  pressure  of  the  flow  {due  to  the  highly  non-isentropic  nature  of  the 
downstream  shocks  encountered}  and  2)  provide  very  accurate  metering  of  the  air  flow 
(when  critical,  flow  is  only  a  function  of  atmospheric  upstream  pressure  and 
temperature).  Design  standards  for  this  type  of  vacuum  “load  chamber”  configuration  are 
available  .  Major  Saupe  can  provide  design  guidance  as  well. 

Although  readily  calculable  (as  enthalpy  is  conserved  for  this  type  of  flow), 
downstream  stagnation  temperature  will  be  very  near  to  that  for  upstream  stagnation 
conditions  (useful  for  first  pass  estimates). 

Flow  is  brought  off  the  end  of  the  ASME  style  critical  flow  nozzle  at  a  90°  angle  to  the 
main  flow  direction.  This  is  to  further  aid  in  breaking  down  the  sonic  (or  possibly 
supersonic)  flow  and  also  aid  in  entropy  production  (further  lower  the  stagnation  pressure 
of  the  flow).  One  key  operational  note  is  that  critical  flow  is  assured  when  these  flow 
nozzles  become  silent  (as  downstream  sound  waves  cannot  propagate  upstream  against 
sonic  flow  in  the  nozzle). 

Vacuum  pumps,  whether  liquid  ring  type,  vane  type  or  ejectors  operate  on  a  specific 
performance  curve.  Lowering  the  mass  flow  rate  to  the  pump  will  typically  result  in 
lower  absolute  suction  pressures.  These  pumps  are  essentially  constant  volume  flow 
machines,  though  volume  flow  capacity  can  change  over  suction  pressure  ranges  due  to 
inefficiencies  within  the  pumps  (typically  at  low  and  high  mass  flow  rates).  OEM 
performance  curves  should  always  be  consulted  when  using  a  vacuum  pump  to  ensure 
performance  will  meet  requirements  over  the  range  of  flows  and  pressures  anticipated 
(see  generic  vacuum  system  curve  shown  in  Eigure  14). 
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Figure  14.  Generic  vacuum  pump  performance  curve.  Typical  industrial  units  are  inches 
of  mercury  absolute  (inHgA)  for  pressure  and  pounds  per  hour  for  flow  rate  (#/hr) 

In  the  configuration  provided  in  Figure  17,  test  cell  pressure  is  provided  by  a  valve 
positioned  downstream  of  the  test  section.  As  the  pump  suction  pressure  is  constant  (the 
result  of  a  constant  mass  flow  rate  fixed  by  the  flow  orifice),  this  valve  will  serve  as  a 
pressure  drop  device  to  adjust  cell  pressure  at  a  desired  level  above  pump  suction 
pressure. 

Although  more  complex,  another  approach  would  be  to  provide  a  recycle  loop  around 
the  pump  (or  even  a  metered  external  air  leak)  such  that  the  pump  sees  more  inlet  mass 
flow  and  will  then  ride  up  its  performance  curve  to  a  new  steady  state  higher  pressure. 
Varying  pump  RPM  is  one  option  as  well  (though  there  may  be  practical  limits  to  this 
approach). 

Due  to  the  very  low  dynamic  head  of  this  system,  actual  flow  velocity  measurements 
using  traditional  head  change  methods  (such  as  pitot  tubes,  venturi  tubes  etc.)  appear 
difficult.  Knowing  the  mass  flow  rate  (as  provided  by  the  flow  nozzles)  will  enable 
velocity  to  be  calculated,  on  an  average  sense,  to  a  fairly  high  degree  of  accuracy.  Local 
velocities  could  possibly  be  measured  by  other  means  such  as  hot  wire  probes  or 
calibrated  turbines.  Mass  flow  rates  in  the  system  can  also  be  cross  checked  against  a 
flow  meter  placed  at  the  vacuum  pump  discharge  and/or  by  using  the  vacuum  pump 
performance  curve. 

4.3.1  Minimum  Vacuum  Pump  Requirements 

For  the  prescribed  flow  conditions  and  geometry,  the  pump  required  will  need  to 
pump  1060  ACFM  (Actual  Cubic  Feet  per  Minute)  of  air  from  73.67  mmHgA  (2.90 
inHgA)  to  760  mmHgA  (29.92  inHgA)  prior  to  accounting  for  line  losses.  Hanscom 
AFRL  currently  has  a  Bush  vacuum  pump  that  can  handle  117  ACFM  at  73.67  mmHgA 
(using  a  7.5  HP  electric  motor).  Initial  estimates  would  require  a  pump  similar  to  that 
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encountered  in  a  municipal  steam  turbine  power  plant  bottoming  cycle,  typically  a  liquid 
ring  single  or  two  stage  pump,  consuming  close  to  150  HP.  Steam  jet  ejectors  could  be 
used  as  well,  and  would  consume  approximately  1700  Ib/hr  of  motive  steam  and  150 
GPM  of  cooling  water 

If  run  times  were  to  be  short,  the  possibility  of  using  a  vacuum  “accumulator”  could 
also  be  investigated.  A  large  vessel  evacuated  to  a  pressure  much  lower  than  76  mmHgA 
could  be  initially  evacuated  over  time  to  provide  a  short  run  of  large  capacity  vacuum 
source. 

Vacuum  capacity  requirements  could  be  marginally  reduced  however  by  operating 
with  the  system  at  a  higher  air  temperature  (i.e.  resulting  in  a  higher  pressure  for  a  given 
density  requirement),  see  Figure  15.  A  chiller  placed  just  upstream  of  the  pump  could 
then  serve  to  lower  the  ACFM  entering  the  pump. 

It  is  anticipated  that  due  to  the  large  costs  associated  with  procuring  and  operating  a 
vacuum  system  of  this  magnitude  this  configuration  may  not  be  desirable. 


Chamber  Pressure  (inHgA) 


Figure  15.  Required  chamber  pressure  vs.  test  cell  temperature  (F) 

4.4  Closed  Loop  Vacuum  Pump  configuration 

The  Closed  Loop  Vacuum  Pump  configuration  is  nearly  identical  to  the  Once  Through 
configuration  with  the  exception  that  instead  of  drawing  from,  and  discharging  to,  open 
atmosphere,  the  working  fluid  is  placed  in  a  closed  loop  configuration  utilizing  the  use  of 
a  retainment  vessel  (see  Figure  18). 

In  this  configuration  it  would  be  possible  to  change  the  working  fluid  to  one  with  a 
much  lower  molecular  weight  (MW)  than  that  of  air.  This  change  would  result  in  a  much 
higher  absolute  test  cell  pressure  and  smaller  vacuum  system  AP  requirements  for  the 
prescribed  density.  Figure  16  illustrates  the  effect  on  pump  suction  pressure  of  lowering 
MW. 
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Figure  16.  Effect  of  molecular  weight  and  temperature  (°F)  on  test  cell  operating 
pressure. 

The  plot  indicates  that  as  MW  decreases,  vacuum  pump  AP  performance  requirements 
are  lowered,  until  for  very  low  MW  (such  as  with  Hydrogen)  chamber  pressure  could 
actually  be  higher  than  atmospheric.  In  such  cases,  a  vacuum  pump  could  be  replaced  by 
a  much  more  economical  centrifugal  pump  or  blower  for  circulation.  Obtaining  and 
retaining  gas  purity  in  such  a  system  may  prove  formidable  however. 

The  same  advantages  remain  in  this  system  with  regard  to  known  flow  rates  via 
critical  flow  nozzles  and  pump  performance  curves.  Instrumentation  to  be  tested  in  the 
cell  should  be  validated  as  suitable  for  whatever  fluid  is  chosen. 


4.5  Closed  Loop  Fan  configuration 

The  Closed  Loop  Fan  configuration  is  shown  in  Figure  19.  In  this  configuration  air  is 
used  as  the  test  fluid  and  is  circulated  via  a  fan  (similar  to  a  closed  loop  wind  tunnel).  Air 
density  is  controlled  through  the  use  of  an  external  vacuum  pump.  In  this  manner  a 
vacuum  pump  is  used  to  initially  evacuate  the  test  system  down  to  a  prescribed  density. 
A  fan  then  circulates  the  flow.  The  vacuum  pump  remains  in  standby  mode  to  remove 
any  air  that  should  it  leak  into  the  system  or  counter  any  increase  in  gas  temperature.  The 
vacuum  pump  used  for  this  application  could  be  rather  small  (provided  air  in-leakage  to 
the  system  was  minimal). 

This  configuration  has  the  advantage  of  low  capital  cost  and  low  cost  of  operation. 
One  disadvantage  stems  from  determining  an  accurate  method  of  flow  measurement. 
The  very  low  dynamic  heads  associated  with  this  type  of  flow  may  require  some  sort  of 
hot  wire  probing  or  low  mass  turbine  insertion  into  the  flow  field. 
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Fan  sizes  and  performance  would  need  to  be  looked  at  closely  in  order  to  ensure 
adequate  performance.  Line  losses  should  also  be  closely  reviewed. 

As  a  first  estimate,  a  circulation  fan  capable  of  handling  1060  ACFM  of  air  would 
typically  have  a  blade  diameter  of  12”.  Capacity  would  be  largely  a  function  of 
backpressure  (line  loss).  If  the  test  section  were  to  be  constructed  of  four  lengths  of 
straight  pipe  and  four  elbows,  line  losses  may  be  on  the  order  of  that  shown  in  Table  2. 

Flow  Conditions: 

-  Air  @  462  #/hr,  70°  F  and  74  mmHgA 
Piping: 

Four  elbows  (L/D  =  20) 

Two  10  foot  sections  of  pipe 
Two  4  foot  sections  of  pipe 

Material: 

-  PVC  piping  and  clear  Lexan 
Dimensions: 

-  12”,  10”  or  8”  I.D.  Main  Piping  w/  8”  I.D.  of  flow  test  section 
Table  2.  Estimate  of  line  losses^ 


Parameter 

12” Main  Pipe 

10” Main  Pipe 

8”  Main  Pipe 

Velocity  in  12”  Pipe  (ft/s) 

22.4 

22.4 

22.4 

Velocity  in  8”  Test  Section  (ft/s) 

50.6 

50.6 

50.6 

AP  per  foot  of  main  pipe  (inH20) 

0.00038 

0.00091 

0.00267 

AP  per  total  pipe  28  ft  (inH20) 

0.01088 

0.02588 

0.07506 

AP  per  Elbow  (inFl20) 

0.00776 

0.01543 

0.03574 

AP  per  4  Elbows  (inH20) 

0.03102 

0.06167 

0.14300 

8”  Test  Section  (est)  (inH20) 

0.02500 

0.02500 

0.02500 

Total  Estimated  Pressure  Drop  (inH20) 

0.06689 

0.11255 

0.24306 

A  typical  12”  circulation  fan,  such  as  Dayton  Fan’s  axial  duct  Model  4TM80  has  the 
following  performance  at  standard  atmosphere  sea  level  conditions: 

1627  ACFM  @  0.000  inH20  backpressure  1182  ACFM  @  0.250  inHiO  backpressure 
1465  ACFM  @  0.125  inH20  backpressure 

As  test  density  is  far  below  that  at  standard  atmosphere  sea  level  conditions,  various 
fan  laws  must  be  employed  to  estimate  capabilities  at  these  lower  densities. 

The  following  “test  condition”  performance  for  this  fan  is  estimated  using  on-line 
conversion  programs  at  http://www.torringtotiresearch.com/fancalc.html 
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Table  3.  Estimate  of  fan  performanee  in  low  density  eonditions 


Parameter 

Sea  Level 
Density 

Test  Density 

Max  allowable  baek  pressure  (mH20)  at  1627  ACEM 

0.00 

0.00 

Max  allowable  baek  pressure  (mH20)  at  1465  ACEM 

0.125 

0.010 

Max  allowable  baek  pressure  (inH20)  at  1182  ACEM 

0.250 

0.020 

Nominal  HP 

0.25 

0.025 

RPM 

1750 

1750 

As  shown  on  Table  3,  it  quiekly  beeomes  apparent  that  mueh  eare  must  be  used  when 
seleeting  a  fan  for  this  duty.  In  this  ease,  a  12”  fan  may  not  be  able  to  overeome  the 
antieipated  pressure  drop  in  the  system.  Alternate  fan  designs,  blowers  or  propellers 
designed  for  high  altitudes  may  need  to  be  employed  in  order  to  provide  the  needed  flow 
rate  at  the  estimated  baek  pressure.  Staging  of  axial  fans  eould  also  be  eonsidered. 
Exeeeding  a  fan’s  rated  maximum  diseharge  pressure  should  be  avoided,  as  stalling  of  the 
blades  eould  be  one  possible  outeome.  One  useful  artiele  on  high  altitude  propeller 
design  ean  be  found  at:  http  ://gltrs.  grc.nasa.gov/cgi-bin/GLTRS/browse.pl?  1 998/TM- 1998- 
206637.html  Artiele  refereneed  is  NASA/TM — 1998-206637. 

Eurther  inereasing  the  main  eonduit  piping  size,  though  eostly,  may  provide  an 
aeeeptable  method  to  lower  system  pressure  drop. 


4.6  Test  Cell  (General) 

The  test  eireuit  should  be  designed  to  withstand  the  atmospherie  forees  against  an 
internal  vaeuum.  Round  piping  is  often  the  most  suitable  eonduit  for  this  type  of  serviee. 
One  material  readily  used  is  thiek  walled  PVC  piping.  Eexan  piping  eould  easily  be  used 
in  the  8”  test  seetion.  These  materials  provide  very  smooth  internal  non-eorrosive 
surfaees,  are  easy  to  work  with  and  are  relatively  light. 

Axial  fans  originally  enelosed  in  OEM  sheet  metal  would  most  likely  need  to  be  re¬ 
installed  into  PVC  piping,  should  the  steel  shell  not  be  rated  for  vaeuum  serviee.  If  non 
axial  blowers  are  employed  a  struetural  assessment  of  the  often  flat  eontainment  surfaees 
would  need  to  be  eondueted. 

Depending  on  the  flow  eonditioning  required,  more  sophistieated  flow  fields  eould  be 
eonstrueted.  One  helpful  wind  tunnel  design  link  is: 

http  ://navier.  Stanford,  edu/bradshaw/tunnel/index.html 
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4.7  Conclusions 


Advantages  and  disadvantages  of  the  three  eonfigurations  investigated  are  provided 
below: 

Once  Through 

Advantages:  Simple  eonstruetion,  easy  to  meter  flow,  low  eost  of  piping. 

Disadvantages:  Highly  ineffieient,  large  vaeuum  souree  required,  large  eost  of  operation. 

Closed  Loop  Vacuum  Pump 

Advantages:  Simple  eonstruetion,  easy  to  meter  flow,  low  eost  of  piping. 

Disadvantages:  Highly  ineffieient,  large  vaeuum  souree  required  (depending  on  fluid), 
high  eost  of  operation,  higher  eost  to  obtain  and  eontain  working  fluid  than  using  air. 

Closed  Loop  Fan 

Advantages:  Simple  eonstruetion,  low  eost  of  eapital  equipment,  low  eost  of  operation. 
Disadvantages:  High  altitude  fan  performanee  needed. 

Due  to  the  above,  the  Closed  Loop  Fan  eonfiguration  appears  to  be  the  most  desirable, 
however,  depending  upon  equipment  availability  and  duration  of  testing,  the  other 
eonfigurations  may  ultimately  prove  more  benefieial. 
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ASME  Critical  Flow  Nozzle(s)  -  to  provide 


Discharge  to 
atmosphere  and 
secondary  flow 
measurement 


Vacuum  Pump  (for 
steady  state,  long  runs) 
or  Vacuum  Source 
(such  as  large  evacuated 
\assel)  for  short  duration 
runs 


Figure  17:  Once  through  configuration 
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Figure  18.  Closed  Loop  Vaeuum  Pump  Configuration 
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Figure  19.  Closed  Loop  Fan  Configuration 
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5.  LOW  DENSITY  TUNNEL  DESIGN  AND  CONSTRUCTION 
5.1  Introduction  and  Background 

Duell  returned  in  the  summer  of  2006,  continuing  the  work  started  in  the  summer  of  2005, 
covered  in  section  4.3,  on  high-sensitivity  flow  measurement  instruments  for  low-density  flows. 
The  next  step  was  to  design  and  build  a  wind  tunnel  to  simulate  high-altitude  conditions.  The 
first  use  of  this  wind  tunnel  would  be  to  investigate  the  frequency  response  of  the  ionic 
anemometer  and  its  sensitivity  to  angle  of  flow. 

Saupe’s  trade  off  study,  covered  in  Section  4,  concluded  that  a  wind  tunnel  based  on  a  closed 
loop  fan  configuration  would  be  the  best  design  for  our  purposes.  An  attempt  to  obtain  an 
estimate  for  time  and  cost  of  a  commercially  manufactured  tunnel  was  made  with  Aerolab; 
however  they  declined  to  bid  on  the  project  citing  a  high  volume  of  projects  over  the  summer 
and  the  uniqueness  of  the  request. 

This  report  contains  four  major  sub-sections;  planned  tasks,  which  includes  the  plan  for  the 
wind  tunnel  and  subsequent  testing;  completed  tasks,  which  details  what  was  completed  this 
summer;  remaining  tasks,  which  covers  what  items  need  to  be  completed  for  tunnel  operation; 
and  finally  future  tasks,  which  contains  several  recommendations  to  improve  the  tunnel  and 
instrument  in  the  future. 


5.2  Planned  Tasks 

The  planned  tasks  for  the  summer  included  design  of  the  wind  tunnel,  analysis  of  the  wind 
tunnel,  acquiring  hardware,  constructing  the  tunnel,  verifying  its  operation,  and  testing  the 
instrument  in  the  tunnel.  The  design  of  the  tunnel  includes  everything  from  the  initial 
determination  of  the  size  of  the  tunnel  to  the  selection  of  fasteners  for  securing  the  fan.  The 
analysis  includes  both  determining  the  back  pressure  in  the  tunnel  and  the  operating  temperature 
of  the  tunnel.  Acquiring  hardware  consists  of  both  purchasing  hardware  needed  to  construct  the 
tunnel,  as  well  as  borrowing  as  much  hardware  as  possible  from  other  personnel  and 
departments.  Construction  of  the  tunnel  involves  the  labor  of  cutting  hardware  to  the  correct  size 
and  assembling  the  tunnel.  Operational  verification  includes  ensuring  that  the  leaks  in  the  tunnel 
are  not  excessive  and  that  it  performs  as  required  for  testing.  Testing  the  instrument  will  begin 
with  tests  to  check  its  angular  and  frequency  response.  During  the  course  of  this  summer,  design, 
analysis,  and  acquisition  were  completed;  construction  was  partially  completed;  verification  and 
testing  have  not  been  started. 
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5.3  Completed  Tasks 

This  chapter  of  the  report  details  what  has  been  eompleted  the  summer  of  2006.  The  major 
eomponents  of  this  ehapter  are  design  deeisions,  engineering  analysis,  eleetronies,  and 
eonstruetion. 


5.3.1  Design  Decisions 

The  test  seetion  diameter  was  sized  to  aeeommodate  the  existing  ionie  anemometer,  whieh  is 
approximately  6  in.  long.  An  8  in.  inner  diameter  pipe  was  seleeted  for  the  test  seetion  to  avoid 
signifieant  boundary  layer  effeets.  With  the  desired  maximum  tunnel  wind  speed  of  30  knots, 
this  ereates  a  required  volumetrie  flow  rate  of  1060  CFM. 

The  fan  seetion  of  the  tunnel  was  sized  based  on  the  minimum  fan  size  to  provide  the  desired 
flow  rate.  The  baek  pressure,  diseussed  below  in  Engineering  Analysis,  proved  to  be  a  diffieult 
speeifieation  to  meet.  After  some  researeh  it  was  found  that  a  12”  diameter  fan  eould  meet  the 
flow  and  baekpressure  requirements. 

In  general  the  design  of  the  tunnel  followed  the  best  praetiees  from  literature.  Following  the 
best  praetiees  from  literature,  a  high  aspeet  ratio  (approximately  6)  honeyeomb  sereen  is  used  in 
the  tunnel  just  before  the  test  seetion  to  reduee  turbulenee.  One  eompromise  to  aid  in 
eonstruetion  is  the  omission  of  sereens;  while  sereens  improve  flow  uniformity,  the  diffieulty  of 
mounting  them  in  the  tunnel  led  to  their  absenee. 

Another  notable  deviation  from  best  praetiees  is  the  design  of  the  diffuser.  Diffusers  generally 
expand  gradually,  often  using  three  legs  of  the  tunnel  for  expansion,  to  minimize  pressure  losses 
and  avoid  separation.  Due  to  the  diffieulty  of  forming  a  gradual  expansion  with  PVC,  and  the 
amount  of  labor  needed  to  form  the  tunnel  from  sheet  metal,  it  was  deeided  to  use  a  short, 
tapered  PVC  eonneetion  for  expansion.  The  relatively  small  expansion  ratio,  2.25,  partially 
mitigates  the  shorteomings  of  rapid  expansion. 

Pipe  flanges  have  been  plaeed  on  both  ends  of  the  fan  and  test  seetions  of  the  tunnel  so  they 
ean  be  removed  and  aeeessed  relatively  easily.  Eaeh  pipe  flange  uses  a  slip  ring  with  four  bolts  to 
seeure  it  in  plaee,  and  a  neoprene  gasket  between  the  flanges  helps  to  reduee  leaks. 

The  instrument  is  mounted  in  the  test  seetion  using  a  eollar  and  set  serews  to  attaeh  the 
instrument  to  a  1”  dia  PVC  tube.  This  tube  passes  through  the  plug  that  seals  off  the  test  seetion, 
and  terminates  with  an  elbow  and  end  eap.  This  allows  the  tunnel  to  be  sealed  off,  while 
allowing  the  instrument  to  be  rotated  for  angular  response  tests.  A  drawing  of  the  mount  is  below 
in  Figure  20. 
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Ionic  Anemometer 

Figure  20.  Drawing  of  Instrument  Mount 

The  vaeuum  pump  is  attaehed  to  the  tunnel  via  1”  PVC  piping.  Two  ball  valves  are  used;  one 
between  the  tunnel  and  the  pump  and  one  between  the  atmosphere  and  the  pump.  The  valve 
between  the  tunnel  and  pump  is  used  to  regulate  the  amount  of  air  being  removed  from  the 
tunnel.  The  seeond  valve  to  the  atmosphere  allows  some  air  into  the  pump  so  it  does  not  run  at 
no-load  for  an  extended  period  of  time.  Several  threaded  eonneetions  have  been  added  to  make 
disassembly,  for  example  for  the  purpose  of  ehanging  the  distanee  from  the  pump  to  the  tunnel, 
easier.  A  diagram  of  the  eonneetion  is  below  in  Figure  21. 
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Figure  21.  Vacuum  Pump  Connection 


A  motor  mount  was  designed  to  hold  the  motor  and  fan  in  the  center  of  the  fan  section  of  the 
wind  tunnel.  The  mount  fits  between  the  motor  and  the  fan  due  to  the  position  of  the  mounting 
holes  on  an  external  rotor  motor.  The  motor  mount  to  hold  the  fan  inside  the  tunnel  is  shown 
below  in  Figure  22. 
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0.25  INCH  ALUMINUM  PLATE 


5.3.2  Engineering  Analysis 

The  first  engineering  analysis  performed  on  the  tunnel  was  to  determine  the  pressure  loss 
around  the  tunnel.  The  worst  test  ease  was  determined  to  be  a  simulation  of  55,000  ft  at  a  speed 
of  30  knots.  The  back  pressure  in  this  case  was  found  to  be  0.25  in-H20  at  test  density.  This  is 
equivalent  to  approximately  2.9  in-H20  at  sea  level  conditions,  which  is  how  fans  are  rated. 
Many  commercial  providers  of  HVAC  equipment  were  contacted  to  attempt  to  find  a  fan  with 
these  performance  characteristics.  The  ELTA  line  of  fans  from  Continental  Fan  was  found  to  be 
suitable  for  this  application.  Following  the  fan  selection,  a  motor  needed  to  be  found  to  power 
the  fan.  The  two  key  requirements  are  that  the  motor  can  vary  its  power  setting  within  a  wide 
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range,  to  create  a  range  of  speeds  within  the  tunnel,  and  the  motor  must  be  brushless,  to  avoid 
potential  problems  with  sparking  at  low  densities.  The  EC90  brushless  flat  motor  from  Maxon 
Motor  was  found  to  be  suitable. 

The  second  engineering  analysis  performed  on  the  tunnel  was  to  determine  the  operating 
temperature.  There  was  some  concern  that  the  temperature  within  the  tunnel  would  become 
excessively  warm  during  prolonged  operation.  However  an  analysis  of  the  heat  flux  through  the 
tunnel  walls  showed  that  the  temperature  inside  the  tunnel  would  be  just  20  F  warmer  than 
outside  the  tunnel  in  the  worst  case  (again,  55,000  ft  simulation  at  30  kts).  The  actual 
temperature  difference  during  operation  is  likely  to  be  even  lower,  since  the  cold  air  leaking  into 
the  tunnel  will  reduce  the  temperature  inside  the  tunnel. 

Heating  of  the  tunnel  air  was  identified  as  a  possible  problem.  Exchange  of  heat  to  the 
ambient  air  through  natural  convection  was  investigated  to  determine  if  any  extraordinary 
measures  would  be  required.  Though  the  thermal  conductivity  of  PVC  is  quite  low,  the  vast 
surface  area  provided  by  the  wind  tunnel  piping  provides  a  substantial  heat  transfer  surface. 

A  radial  heat  transfer  model  was  generated  using  forced  convection  on  the  internal  surfaces, 
conduction  through  the  wall  of  the  piping,  and  then  natural  convection  to  the  surrounding  air. 
Internal  forced  convection  was  modeled  using  the  classic  Dittus-Boelter^"^  formulation: 

Nu  =  0.023  Re/®  Pr“'  (24) 

Where  Nu  is  the  Nusselt  number,  Nu  =  hD/K  ,  h  is  the  heat  transfer  coefficient,  D  is  diameter, 
and  K  is  the  thermal  conductivity.  Pr  is  the  Prandtl  Number,  Vx  =  vl a  ,  v  =  kinematic  viscosity 
and  a  =  thermal  diffusivity.  Re  is  the  Reynolds  Number,  Re  =  ED/v  ,  where  V  is  velocity. 
Conduction  within  the  pipe  sections  took  the  standard  cylindrical  conduction  form: 

(25) 


Qcond  — 


IttKL 


Natural  convection  was  modeled  using  the  following  relation  for  externally  buoyant  convection 
on  a  horizontal  cylinder 


Nu  =  ^ 


0.6  + 


0.38(GrPr)/  l  +  (0.559/Pry 


16 


16/ 


(26) 


Where 


Gr  = 


D^gaAT 


(27) 


With  a  known  heat  flux  through  each  section  of  pipe  and  a  dependence  upon  the  temperature 
difference,  AT,  for  the  ambient  Nu  number,  an  iterative  scheme  was  used  to  arrive  at  a  complete 
solution. 
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The  results  of  this  investigation  for  the  heaviest  power/thermal  loading  (55  000  ft),  show  a  4 
degree  rise  from  ambient  temperature  to  the  outer  pipe  wall  surfaee,  a  nominal  inerease  of  6 
degrees  aeross  the  pipe  material  itself  and  another  8  degree  rise  from  the  pipe  wall  inside  surfaee 
temperature  to  the  bulk  fluid  flow  temperature.  This  is  approximately  an  18  degree  rise  above 
ambient  (or  an  approximate  bulk  flowing  temperature  of  90°F  based  upon  ambient  eooling). 

If  the  above  ealeulations  are  reasonably  aeeurate,  no  aetive  eooling  sehemes  (i.e.  foreed  air 
leaks)  should  be  required  (though  they  eould  be  readily  employed  if  needed).  Room  fans  eould 
also  be  utilized  to  assist  with  ambient  eooling  if  needed. 

5.3.3  Electronics 

To  provide  usable,  reeordable,  signal  levels,  eurrent  flow  indueed  in  the  eolleetion  rods  of  the 
anemometer  must  be  proeessed  utilizing  analog  signal  eonditioning  eleetronies.  A  simple 
sehematie  of  the  eireuit  is  shown  in  Figure  23.  Current  from  eaeh  eolleetor  is  eonverted  to 
voltage,  and  the  sum  and  differenee  of  the  voltages  are  eomputed  utilizing  operational  amplifier 
based  eireuitry.  Ineluded  in  eaeh  of  these  eireuits  is  an  offset  adjustment  to  eompensate  for 
geometrie  ineonsisteneies  with  probe  eonstruetion,  as  well  as  any  eleetronie  input  offsets 
inherent  with  operational  amplifier  based  eireuits. 

The  sum  and  differenee  signals  are  then  amplified  to  provide  appropriate  voltage  levels 
allowing  maximum  resolution  from  a  PC  based  data  aequisition  system.  Due  to  the 
suseeptibility  of  the  anemometers’  open  eolleetors’  to  ambient  eleetronie  noise;  it  is  also 
neeessary  to  provide  filtering  to  mitigate  signals  indueed  by  a  laboratory  setting  (60  Hz  noise). 
In-flight  operation  away  from  any  AC  powered  equipment  may  not  require  this  step.  These 
amplified  and  filtered  signals  are  then  reeorded  utilizing  a  16-bit  PC  based  data  aequisition 
system  eapable  of  .3  mV  resolution  over  a  +/-  lOV  input  range. 


Figure  23.  Simple  sehematie  of  filtered  sum  and  differenee  eireuit.  Current  from  eaeh  eolleetor 
eomes  in  from  the  left.  The  filtered  sum  is  the  upper  output,  the  filtered  differenee  on  the 
bottom. 
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Lt  into  sections  ana  me  iittmgs  were  an  attac 
i-built”  drawing  of  the  wind  tunnel  layout  is 


6.  TUNNEL  INITIAL  INSTRUMENTATION  CHECK  AND  TUNNEL 
CHARACTERIZATION 


In  December  2006,  the  wind  tunnel  was  commissioned  for  sea-level  use.  The  first  part  of  the 
commissioning  process  involved  assessment  and  curve  fitting  of  calibration  information  of 
pressure  and  temperature  instrumentation. 


The  two  primary  pressure  instruments  utilized  were  absolute  in  nature  with  an  applicable 
range  from  standard  atmospheric  pressure  (760  mmHgA)  down  to  full  vacuum  (0  mmHgA). 
These  instruments  were  checked  at  near  full  vacuum  (0.01  mmHgA)  against  an  AFRL  laboratory 
vacuum  measurement  standard.  The  instruments  were  then  checked  against  local  atmospheric 
pressure  measurement  (755  mmHgA),  as  measured  from  an  AFRL  barometric  pressure  standard. 
As  a  cross  check,  the  local  atmospheric  pressure  was  compared  to  the  local  weather  agency 
reported  atmospheric  pressure  (corrected  for  lab  room  elevation  and  atmospheric  temperature). 
Both  methods  produced  excellent  agreement.  As  the  pressure  instruments  used  on  the  tunnel  are 
designed  with  a  linear  response,  the  full  vacuum  and  atmospheric  pressure  points  were  used  to 
create  a  two  point  linear  calibration  curve  for  each  device  as  shown  in  Figure  25. 


**  BLUE  CELLS  ARE  INPUTS 

Pressure  Correction  Cross  Check 


29.81  inHgA  = 
760.00  mmHgA  = 


757.21 

29.92 


mmHgA 

inHgA 


Reported  Barometric  Pressure  = 

757.21 

mmHgA 

Reported  Altimeter  Setting  = 

29.97 

InHgA 

Inst  Height  Above  Sea  Level  = 

275 

ft 

Inst  Height  Above  Sea  Level  = 

275 

ft 

Atmospheric  Temp  = 

45 

°F 

Pressure  Correction  = 

0.294 

inHgA 

Exponent  = 

0.0102195  mmHgA 

Corrected  Local  Pressure  = 

29.676 

inHgA 

STD  DAY  CORRECTION  (Aviation) 

Corrected  Local  Pressure  = 

749.51 

mmHgA 

Corrected  Local  Pressure  = 

753.81 

mmHgA 

Pressure  Elevation 

A  Correction 

29.00  863 

0.92 

Note:  Tsat  @  30  mmHgA  =  82  Deg  F 

29.01  853 

0.91 

29.02  844 

0.90 

29.03  834 

0.89 

Serial  Number 

SS13097/FA160 

Serial  Number 

TT10557  /  FA160 

Reading 

Real 

Fit 

Reading 

Real 

Fit 

mmHgA 

mmHgA 

mmHgA 

Reading 

Corrected 

mmHgA 

mmHgA 

mmHgA 

760 

755 

755.0 

mmHgA 

mmHgA 

800 

755 

755.0 

28 

0.01 

0.01 

758 

752.93 

30 

0.01 

0.01 

Correction  Curve: 

Correction  Curve: 

Real  =  (Readinq)*1.0314  -  28.869 

Real  =  (Readinq)*0.9805  -  29.405 

Pressure  Correction  Curve 

y=0.9805x- 29.405 
=  1 


Figure  25.  Pressure  gage  calibration  data  and  curve  fits 
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OMEGA  Thermistor 


Resistance:  10,000  ohms  @°25 
10to8tfC  =  50to176F 


[Conversion  Equation 

=(-5.83475106E-19)*Ohms'^5  +  (3.36817467E-14)*OhmsM  -  (7.54527227E-1Q)*Ohms^3 


+  (8.40899264E-06)*Ohms''2  -  (5.220441 91  E-02)*Ohms  +  2.34641272E+02 
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Figure  26.  Thermistor  ealibration 

Tunnel  temperature  was  determined  using  a  thermistor.  This  thermistor  was  provided  with  a 
ealibration  eertificate  from  the  manufacturer,  relating  temperature  (°C)  to  resistance  (ohms) 
through  the  use  of  tabulated  data.  As  the  response  was  highly  non-linear,  a  5**^  order  polynomial 
was  fit  to  the  data.  The  details  are  provided  in  Figure  6. 

Initial  runs  of  the  tunnel  focused  on  determining  the  fan  response  to  tunnel  back  pressure 
(effectively  trying  to  calibrate  the  fan  RPM  to  a  measured  mean  tunnel  velocity).  During  very 
low  density  conditions,  our  velocity  check  device,  a  pitot  tube,  would  be  of  very  little  use  to 
evaluate  tunnel  speed  (as  our  A  pressure  instrumentation  did  not  have  the  accuracy  required  for 
such  low  AP  operating  conditions).  Applying  fan  laws  (for  density  changes)  and  expected  fan 
curve  shifts  (for  AP)  would  be  a  basis  for  determining  tunnel  velocity  at  high  altitude  conditions. 
Mean  tunnel  velocity  as  a  function  of  fan  RPM  and  pressure  drop  is  shown  in  Figure  27. 
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Figure  27.  3500  RPM  Sea  Level  Fan  Performanee 


Initial  test  run  data  (provided  in  Table  4),  indicated  a  large  discrepancy  between  calculated 
velocity  from  fan  RPM  and  measured  values  from  the  pitot  system  45%  error).  This  error  was 
later  found  to  be  attributed  to  an  error  in  calculating  fan  RPM  as  a  function  of  reported  fan  motor 
frequency  (Hz).  Tests  were  re-run  to  calibrate  the  tunnel  at  sea-level  using  a  high  accuracy  wind 
turbine  meter.  A  velocity  vs  fan  RPM  calibration  curve  is  shown  in  Figure  28. 


Data  were  later  collected  at  138  millibars  (highest  attainable  steady  state  altitude)  using  fan 
frequency  as  the  prime  indicator  of  mean  tunnel  velocity.  As  accurate  pressure  drop  data  was 
not  attained,  it  has  been  assumed  that  the  Velocity  vs  fan  RPM  calibration  curve  for  sea-level 
remained  valid  at  138  millibars. 
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Table  4;  Input  for  Initial  Characterization  Run 


Date:  18-Dec-06 


TP1 

TP2 

TP3 

TP4 

TP5 

TP6 

TP7 

TP8 

TP9 

TP10 

TP11 

TP12 

TP13 

Input  Values 

Time 

12:20 

12:25 

12:30 

12:35 

12:36 

13:36 

14:36 

15:36 

16:36 

17:36 

18:36 

19:36 

20:36 

Trigger  Voltage  (V) 

0.02 

0.50 

0.98 

1.25 

1.50 

1.95 

2.10 

2.40 

2.70 

3.00 

3.10 

3.52 

4.01 

Fan  Hz 

16.00 

49.00 

88.00 

110.00 

130.00 

165.00 

177.00 

202.00 

225.00 

250.00 

258.00 

291.00 

329.00 

Driving  Voltage  (V) 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

50.00 

Driving  Amps 

0.00 

0.01 

0.01 

0.25 

0.20 

0.30 

0.30 

0.40 

0.50 

0.60 

0.70 

0.90 

1.25 

Driving  Power  (Watts) 

0.05 

0.50 

0.50 

12.50 

10.00 

15.00 

15.00 

20.00 

25.00 

30.00 

35.00 

45.00 

62.50 

Fan  Speed  (RPM) 

80 

245 

440 

550 

650 

825 

885 

1,010 

1,125 

1,250 

1,290 

1,455 

1,645 

RPM  Based  Upon  Fan  Curve  Power 

204 

439 

439 

1,284 

1,192 

1,365 

1,365 

1,502 

1,618 

1,719 

1,810 

1,968 

2,196 

aP  Across  Fan  (inH20) 

0.00 

0.01 

0.03 

0.05 

0.06 

0.11 

0.12 

0.15 

0.19 

0.24 

0.26 

0.25 

0.35 

aP  Across  Pitot  (inH20) 

0.00 

0.00 

0.01 

0.02 

0.02 

0.03 

0.04 

0.05 

0.06 

0.08 

0.08 

0.10 

0.12 

Resistance  Across  Thermistor  (Ohms) 

11,400 

11,440 

11,400 

11,400 

11,400 

11,400 

11,400 

11,400 

11,400 

1 1 ,400 

11,400 

11,400 

1 1 ,400 

Tunnel  Absolute  Pressure  (mmHgA) 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

762.0 

Output  Values 

Fan  Curve  Volumetric  Flow  (ft^/min) 

46 

140 

251 

312 

368 

463 

495 

562 

622 

685 

705 

796 

886 

Pitot  Tube  Volumetric  Flow  (ft^/min) 

0 

70 

140 

177 

207 

258 

279 

312 

348 

388 

395 

442 

488 

Corrected  Tunnel  Pressure  (mmHgA) 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

757.06 

Tunnel  Temperature  (°F) 

71.01 

70.83 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

71.01 

Density  (Ib/ft^) 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

0.0745 

8"  Flow  Area  (ft^) 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

0.347 

Tunnel  Velocity  (ft/s)  Fan 

2.20 

6.72 

12.02 

14.98 

17.64 

22.21 

23.77 

26.96 

29.84 

32.86 

33.81 

38.20 

42.53 

Tunnel  Velocity  (ft/s)  Pitot 

0.00 

3.35 

6.70 

8.48 

9.94 

12.35 

13.40 

14.98 

16.68 

18.59 

18.95 

21.19 

23.40 

%  Deviation 

N/A 

50.15 

44.26 

43.41 

43.66 

44.37 

43.63 

44.43 

44.10 

43.43 

43.95 

44.53 

44.97 

Viscosity  (Ib-s/ft^) 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

3.83E-07 

8"  Section  Reynolds  #  (Fan  Curve) 

8.83E+03 

2.70E+04 

4.83E+04 

6.02E+04 

7.09E+04 

8.92E+04 

9.55E+04 

1.08E+05 

1.20E+05 

1.32E+05 

1.36E+05 

1.53E+05 

1.71E+05 

8"  Section  Reynolds  #  (Pitot) 

O.OOE+00 

1.35E+04 

2.69E+04 

3.40E+04 

3.99E+04 

4.96E+04 

5.38E+04 

6.02E+04 

6.70E+04 

7.47E+04 

7.61  E+04 

8.51E+04 

9.40E+04 

Tunnel  Friction  Factor  (K)  (Fan  Curve) 

N/A 

0.00017719 

0.000194 

0.000201 

0.000201 

0.000215 

0.000209 

0.000211 

0.000211 

0.000223 

0.000227 

0.000171 

0.000191 

Tunnel  Friction  Factor  (K)  (Pitot) 

N/A 

0.00071307 

0.000624 

0.000627 

0.000633 

0.000694 

0.000657 

0.000682 

0.000675 

0.000697 

0.000724 

0.000557 

0.00063 

Tunnel  Watts  (Fan  Curve  Velocity) 

0.000 

0.132 

0.824 

1.650 

2.699 

5.762 

6.866 

10.097 

13.735 

19.389 

21.519 

23.377 

35.918 

Tunnel  Watts  (Pitot  Velocity) 

0.000 

0.066 

0.459 

0.934 

1.520 

3.206 

3.871 

5.611 

7.678 

10.969 

12.062 

12.967 

19.765 

8.0 

7.0 

V 

e  6.0 

p  4.0 
t  3.0 

(m/s)  2  ° 
1.0 
0.0 


y  =  0.0206439X 
R2=  0.9999762 


0 


50 


100  150  200  250  300 

Fan  Motor  Frequency  (Hz) 


350  400 


Figure  28.  Sea  level  fan  calibration  and  curve  fit 
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7.  IONIC  ANEMOMETER  ONE-ATMOSPHERE  WIND  TUNNEL 
TESTING 

The  results  for  our  sea-level  wind  tunnel  runs  are  shown  below  in  Figures  29  and  30.  As  can 
be  seen  in  Figure  29,  the  ratio  of  the  difference  current  to  the  sum  current  of  the  ionic 
anemometer  increases  linearly  with  windspeed.  This  is  as  expected  from  Equation  23.  Figure  30 
shows  the  comparison  of  the  calculated  windspeed  based  on  the  data  collected  by  the  ionic 
anemometer,  and  the  actual  measured  windspeed  of  the  tunnel  using  a  conventional  propeller 
anemometer.  However,  in  conducting  these  tests,  the  method  in  which  the  ionic  anemometer 
data  was  collected  does  cause  a  slight  amount  of  uncertainty.  To  take  a  reading  at  a  given 
windspeed,  first  the  variable  input  voltage  for  the  fan  motor  was  adjusted  to  the  desired  amount. 
Afterwards,  the  windspeed  was  allowed  to  propagate  in  the  wind  tunnel  for  several  seconds  until 
a  steady  state  was  achieved.  Despite  this,  there  were  still  small  fluctuations  in  windspeed  that 


Windspeed  (kts) 


Figure  29.  Windspeed  deduced  from  fan  speed  vs  the  ratio  of  collector  current  difference  to 
collector  current  sum  tEa.  23')  for  1  atmosnhere  nressure  in  the  wind  tunnel. 


Windspeed  (kts) 


Eigure  30.  Comparison  of  windspeed  deduced  from  the  calibrated  ionic  anemometer 
(blue  circles)  to  the  speed  measured  by  a  propeller  anemometer  (pink  diamonds),  as  a 
function  of  the  A/Z  current  ratio. 
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the  very  sensitive  ionie  anemometer  was  able  to  pick  up.  The  recorded  values  for  the  ratios  were 
all  average  values  over  several  seconds  while  at  steady  state.  This  will  become  more  significant 
during  the  higher  altitude  simulations. 

8.  LOW  DENSITY  TESTING 

The  main  problem  we  encountered  with  the  wind  tunnel  was  maintaining  a  steady  pressure  for 
the  desired  altitude  we  were  simulating.  The  original  plan  was  to  take  data  at  a  pressure  of  90 
millibars.  However,  this  introduced  too  much  uncertainty  in  our  data  collection.  This 
uncertainty  was  primarily  caused  by  small  leaks  in  our  system.  As  mentioned  previously,  the 
activation  voltage  of  the  ionic  anemometer  is  dependant  upon  the  tunnel  density.  At  much  lower 
pressures,  the  activation  voltage  of  the  anemometer  is  also  much  less.  The  leaks  in  the  tunnel 
preventing  the  pressure  from  maintaining  90  millibars  were  affecting  the  readings  to  the  point 
where  we  were  collecting  not  only  wind  velocity  perturbations,  but  also  variations  due  to 
activation  voltage  variances.  Ultimately,  we  decided  to  collect  the  high  altitude  data  at  138 
millibars.  This  allowed  us  to  simulate  high  altitude  performance  while  keeping  the  pressure  in 
the  tunnel  steady  enough  to  avoid  activation  voltage  uncertainty. 

Once  it  was  decided  to  simulate  the  high  altitudes  at  138  millibars,  the  same  windspeed  versus 
ratio  test  was  run.  The  result  of  this  is  shown  in  Figure  3 1 .  While  the  ratio  is  still  increasing 
with  windspeed,  it  is  not  nearly  so  linear  as  the  one  atmosphere  test.  Even  though  tunnel 
pressure  was  more  stable  than  at  90  millibars,  there  were  enough  leaks  to  make  averaging  our 
ratio  values  over  a  test  run  difficult.  We  collected  all  data  within  the  range  of  experienced  ratios, 
but  the  fact  that  our  activation  voltage  was  still  affecting  our  readings  means  the  error  bars  are 
much  greater  than  what  would  be  encountered  at  sea  level.  Despite  this,  the  fact  that  the  ionic 
anemometer  was  getting  increased  ratios  in  such  a  way  with  increasing  windspeed  is  very 


Windspeed  (kts) 


Figure  31.  Windspeed  based  on  fan  speed  vs.  A/Z  current  ratio  at  138  millibars 
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promising  for  the  argument  of  using  it  at  high  altitudes,  espeeially  sinee  pressures  would  be 
nearly  eonstant  for  an  airship  at  constant  altitude. 

The  last  test  was  perhaps  the  most  important,  actual  turbulence  measurement.  Data  was 
recorded  for  8  second  periods.  The  conditions  at  the  time  of  the  data  collection  were:  P  =  138 
mbars,  T  =  22°C.  This  is  a  density  of  0.163kg/m^,  which  is  approximately  16.2  km  or  53  000  ft. 
Due  to  the  presence  of  turbulent  flow  in  the  tunnel,  this  test  not  only  checked  the  response  time 
of  the  ionic  anemometer  (a  major  limiting  factor  of  the  propeller  anemometer),  but  also  provided 
data  to  perform  statistical  analysis.  When  the  ratios  were  collected,  they  were  converted  to  the 
corresponding  windspeed  in  knots  as  shown  in  Figure  32,  a  typical  time  series.  The  figure  shows 
over  4000  points  of  data  collected  in  8  seconds  of  testing,  with  speeds  varying  from 
approximately  7.9  to  8.7  knots.  Clearly,  the  ionic  anemometer  has  a  very  rapid  response. 


Data  ID 


Figure  32.  Time  series  of  velocity  data  from  ionic  anemometer.  The  x  axis  is  the  number  of  the 
sample,  the  vertical  axis  is  the  speed  in  knots. 
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9.  IONIC  ANEMOMETER  TURBULENCE  ANALYSIS 


The  data  from  the  ionic  anemometer  test  is  essentially  a  time  series  of  velocity  measurements, 
which  can  be  analyzed  for  tunnel  turbulence  levels.  The  time  series  is  shown  in  Figure  32.  The 
objective  of  the  analysis  was  to  verify  that  the  turbulence  levels  were  reasonable.  Ideally,  one 
would  have  velocity  measurements  from  a  calibrated  source  to  compare  to  the  anemometer 
output.  This  was  not  the  case.  A  hot  wire  probe  would  be  one  such  source  of  data,  but  a  probe 
that  could  be  installed  through  the  vacuum  wall  was  not  available  in  the  time  of  the  run.  Instead 
the  anemometer  data  were  subjected  to  statistical  analysis  and  tested  against  the  reasonableness 
of  the  results. 

The  initial  calculations  were  for  mean  velocity  and  the  root-mean-square  (rms)  of  the  velocity 
fluctuations.  The  results  of  the  analysis  were  that  the  mean  tunnel  velocity  was  4.24  m/s  (8.25 
knots),  and  the  rms  of  the  variance  (or  standard  deviation  of  the  measurement)  was  0.0738  m/s 
(0.1435  knots).  The  ratio  of  rms  to  mean  velocity  gives  the  turbulence  intensity,  which  is  0.0174 
or  1.74%.  This  is  considered  to  be  on  the  lower  end  of  typical  turbulence  for  large  pipes  and 
ventilation  flows,  which  are  typically  between  1%  and  5%. 

Next  the  results  were  analyzed  using  power  spectral  density  (PSD)  and  structure  functions. 
The  PSD  is  defined  as  the  Fourier  Transform  of  the  data  times  its  conjugate^^.  The  PSD  of  the 
velocity  data,  created  using  a  fast  Fourier  transform  computer  algorithm,  is  shown  in  Figure  33. 
The  red  line  is  the  slope  hypothesized  by  Kolmogorov^’  for  the  PSD  of  homogeneous  isotropic 
turbulence  for  the  scales  between  the  outer  scale  on  the  low  frequency  side  to  the  inner  scale  on 
the  high  frequency  side.  The  outer  scale  is  the  size  of  the  largest  turbulent  eddies  that  form  from 
a  disturbance,  and  the  energy  of  turbulent  motion  is  postulated  to  be  conserved  as  it  cascades 
down  to  the  inner  scale,  where  the  energy  begins  to  be  transformed  to  heat  through  viscous 
processes.  Also  shown  is  the  nominal  frequency  of  a  low  pass  filter  used  to  minimize  60Hz 
interference  known  to  be  present  in  the  laboratory.  Data  to  the  right  of  the  filter  line  should 
show  attenuation,  until  the  noise  floor  is  reached,  where  the  PSD  should  level  out.  The 
transformation  from  the  frequency  domain  to  the  spatial  domain  is  accomplished  using  the  mean 
flow  rate  of  the  tunnel,  and  assuming  “frozen  flow”,  i.e.  the  eddies  are  assumed  to  not  change  as 
they  flow  past  the  sensor.  Typically,  the  PSD  for  turbulence  is  plotted  as  a  function  of  spatial 
frequency  k  (  m'^),  k  =  f  iKjV ,  or  about  a  factor  of  1 .5  times  the  temporal  frequency  values  in 
Figure  33  for  our  mean  velocity  of  4.24  m/s. 
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Figure  33.  Power  Speetral  Density  (PSD)  of  the  tunnel  veloeity  data  from  the  Ionic 

Anemometer. 

With  all  the  variation  in  this  PSD  about  all  we  can  say  here  is  that  there  appears  to  be  a  region  of 
Kolmogorov  Turbulence  in  the  data.  The  PSD  does  seem  to  show  a  decrease  in  energy  to  the 
right  of  the  filter  line,  and  suggests  that  the  data  become  dominated  by  noise,  where  the  PSD 
begins  to  level  out  above  about  lOOHz.  This  is,  in  fact  right  about  where  one  would  expect  to 
lose  information,  with  a  1cm  gap  between  the  electrodes  at  4.25  m/s.  The  right  side  of  the  PSD 
ends  at  the  Nyquist  Frequency,  which  is  the  highest  frequency  available  given  the  data 
acquisition  rate,  requiring  that  there  is  at  least  two  data  points  per  cycle.  The  outer  scale  is 
probably  12  to  8  inches,  the  pipe  diameter  before  and  after  the  constriction.  This  translates  to  a 
frequency  of  20  to  14  Flz,  but  it  does  not  appear  to  be  an  obvious  place  where  the  PSD  diverges 
from  the  Kolmogorov  slope. 

Structure  functions  are  another  common  technique  used  to  examine  turbulence  data  .  The 
mathematical  definition  of  a  velocity  structure  function  is: 

(2*) 

Where  (  )  denotes  an  average  over  a  path  s  of  velocity  at  two  points  which  are  a  distance  r  apart. 
The  Kolmogorov  hypothesis  when  expressed  as  a  structure  function  becomes: 

=  (29) 

where  C^j  is  the  velocity  structure  constant.  The  structure  function  of  the  velocity  data  is  shown 
in  Figure  34.  Note  the  smoothing  effect  of  the  structure  function  process. 

On  the  larger  lengths  of  the  structure  function,  on  the  right  side,  it  is  clear  where  the  velocity 
diverges  from  the  Kolmogorov  slope  when  the  outer  scale  is  reached.  In  this  case,  the  outer 
scale  size  agrees  with  the  expected  order  of  magnitude  of  the  tunnel  diameter,  which  reduces 
from  .3  m  (12  in)  to  0.2  m  (8  in)  immediately  before  the  instrument  test  section. 
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The  inner  seale,  where  the  shear  in  the  eddies  is  strong  enough  to  begin  dissipation  into 
thermal  energy,  depends  on  the  strength  of  the  turbulence  and  on  the  kinematic  viscosity  of  the 
fluid From  this  reference,  an  approximate  formula  for  the  inner  scale  is 

Where  v  is  the  kinematic  viscosity  and  s  is  the  kinetic  energy  dissipation  rate,  computed  using 

(31) 

where  a  for  this  equation  is  a  constant  approximately  equal  to  2.  is  determined  by 
extrapolating  the  Kolmogorov  segment  to  the  Im  location  (i.e.  10°),  where  it  is  approximately 
0.025m'^'^^^Vs^.  Next,  using  Equation  31,  s  is  estimated  as  0.0013  m^/s^.  The  kinematic 
viscosity  depends  weakly  on  fluid  temperature  and  strongly  on  density.  For  our  conditions,  the 
viscosity  is  approximately  8.8x10'  m  /s.,  which  gives  an  inner  scale  of  approximately  0.02  m. 


Wind  Tunnel  Turbulence:  Structure  Function 


Figure  34.  Structure  function  of  the  tunnel  velocity  data  from  the  Ionic  Anemometer. 

With  the  structure  function  manifestation  of  the  time  series,  it  does  appear  that  there  is  a 
region  of  turbulence  that  follows  Kolmogorov’s  hypothesis  from  the  30Hz  filter  to 
approximately  the  size  of  the  pipe  diameter.  FInfortunately,  the  30Flz  filter  blocks  the  region  of 
the  structure  function  that  would  have  shown  us  the  inner  scale,  and  higher  priority  laboratory 
activities  prevented  further  effort  on  this  experiment.  With  the  accumulated  information  from 
both  PSD  and  structure  function  plots,  these  analyses  give  us  a  reasonable  confidence  that  the 
ionic  anemometer  is  giving  accurate  velocity  fluctuation  information. 
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10.  FLIGHT  PAYLOAD  CONCEPT 


Upon  completion  of  instrumentation  selection  for  high  altitude  operation,  packaging  and 
deployment  of  the  atmospheric  sensor  package  (HAA  payload)  was  of  prime  concern.  Figure  36 
provides  a  sketch  of  the  concept.  The  purpose  of  the  payload  is  to  measure  mean  temperature, 
pressure,  wind  speed  and  direetion,  instrument  position,  veloeity,  and  acceleration  at  both  the 
HAA  altitude  and  a  position  of  at  least  300  m  below  the  airship.  The  payload  will  weigh 
approximately  10  lbs  with  a  power  of  30  W  at  15  V,  will  have  a  data  uplink,  and  have  the  ability 
to  be  deployed  and  retrieved. 


10.1.  Requirements 

The  payload  package  has  many  design  requirements  associated  with  the  atmospheric 
conditions  at  airship  altitudes  -  not  covered  here.  Here,  we  mention  those  related  to  data 
eollection  and  physieal  design.  For  the  data  collection,  the  driving  force  of  the  design  is  that  it 
must  be  able  to  measure  three  axis  mean  and  fluctuating  wind  velocities.  This  will  primarily  be 
the  ionic  anemometer  (fluctuating),  and  propeller  anemometer  (mean),  and  their  respective 
support  eomponents.  The  package  also  must  measure  mean  and  fluetuating  pressure  and 
temperature.  Finally,  it  must  have  the  ability  to  resolve  relative  motion  between  the  payload  and 
the  host.  For  the  physical  requirements,  the  first  thing  to  consider  is  storage  and  deployment 
from  the  host.  Ideally,  the  payload  would  have  an  attaching  mechanism  to  allow  it  to  be 
deployed  and  retracted  directly  from  the  HAA  and  lowered  into  position  with  an  expanding 


Figure  35.  Payload  eoneept  for  atmospheric  turbulence  sensor  payload  suspended  below  the 
High  Altitude  Airship. 
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(unfolding)  structure  during  deployment  using  a  winching  system.  The  payload  would  have  to 
be  deployed  outside  of  the  HAA  air  stream  influence  (approximately  300  m). 

10.2.  Data  Storage  and  Transfer 

The  payload  will  be  on  station  for  potentially  long  durations  (weeks),  so  it  will  be  important 
to  have  data  storage  and  transfer  capabilities.  Ultimately,  the  data  collected  by  the  payload  will 
be  relayed  to  the  HAA  for  initial  storage  and  transmission  to  the  user,  however,  this  interface  can 
only  be  determined  once  details  of  a  final  host  system  are  available  and  reviewed.  Figure  36 
shows  this  relationship. 

10.3.  Design 

We  propose  that  the  payload  be  lowered  from  the  HAA  via  a  thin  power/data  cable.  Attaehed 
to  this  cable  would  be  an  interface  between  the  harness  that  holds  the  payload  and  the  suspended 
power  cable.  Inside  this  connection  mechanism  would  be  housed  a  data  repeater,  power  supply, 
media  storage  and/or  IR  data  interface.  To  minimize  drag  and  vibration  on  the  suspended  cable, 
a  eontinuous  wrap  of  wire  may  be  necessary  to  break  up  shedding  vortices  should  they  oecur. 
This  wire  would  be  similar  to  the  design  eurrently  being  used  in  late  model  car  antennas.  A  four 
line  harness  would  be  used  to  connect  the  platform,  to  significantly  aid  the  stability  on  the  pitch 
and  roll  axis.  The  most  plausible  material  to  use  for  this  harness  would  be  150  lb  Spectra  line, 
which  is  typical  in  the  performance  kite  arena,  and  of  which  drag  would  be  minimal. 


The  payload  structure  resembles  an  upside  down  glide,  without  an  elevator.  A  vertical 
stabilizer  slab  would  be  located  on  the  lower  surface  empennage  to  allot  the  payload  to  maintain 
a  direetional  reference  with  the  wind.  The  wings  would  have  a  negative  angle  of  attack, 
increasing  the  apparent  weight,  to  help  align  the  payload  under  the  HAA  rather  than  being  blown 
further  downstream.  The  ‘fuselage’  is  circular  in  cross  section  and  can  be  used  to  house 
instrumentation  interfaces,  instrumentation  power  supply,  data  storage  and/or  an  IR  data 
transmitter  and  receiver.  The  payload  ean  be  designed  in  a  manner  that  would  allow  it  to  be 
folded  and  then  auto  erected  when  lowered  from  the  HAA.  Front  and  side  views  of  the  payload 
are  shown  in  Figures  37  and  38. 


Wind 

Velocities 

- ► 

- ► 

- ► 


Payload 

Figure  36.  Sehematic  of  Data  and  Power  paths  for  the  atmospherie  payload  of  the  HAA  or 
any  stratospheric  airship. 
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Figure  37.  Front  view  of  payload  figure  38.  Side  view  of  payload 


10.4.  Instrumentation 

As  mentioned  in  the  requirements  for  the  payload,  there  are  many  different  types  of  data  that 
must  be  eolleeted  with  this  payload.  The  first,  and  most  important,  is  wind  veloeity.  Fluctuating 
velocities  will  be  collected  from  at  least  two  4-collector  ionic  anemometers  oriented  in  planes 
normal  to  each  other,  with  redundancy  in  the  upstream  axis.  Mean  velocities  will  be  collected 
from  both  the  ionic  anemometers  and  a  propeller  anemometer,  which  will  serve  as  real-time 
velocity  calibration.  Temperature  will  be  measured  from  high  accuracy  resistance  temperature 
detectors  (RTDs),  thermistors,  and/or  thermocouples.  The  pressure  value  will  be  obtained  from 
the  HAA’s  onboard  pressure  sensors,  or  included  in  the  payload,  space  and  weight  permitting. 
Finally,  absolute  position  and  position  relative  to  the  FIAA  will  be  obtained  utilizing  Global 
Positioning  System  (GPS)  and  Inertial  Navigation  System  (INS). 
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11.  CONCLUSIONS  AND  RECOMMENDATIONS 


Measuring  velocity  fluctuations  on  a  stratospheric  airship  such  as  the  High  Altitude  Airship, 
at  high  altitude  and  low  speed,  is  a  difficult  task  for  instruments  designed  to  perform  that 
function  at  sea  level  and  low  altitude  flight.  While  variants  of  the  pitot-static  probe  have  been 
used  very  successfully  at  lower  altitudes,  the  combination  of  low  density  and  low  speed  results  in 
dynamic  pressure  differences  that  are  barely  detectible,  especially  for  long  endurance  operational 
instruments.  The  propeller  anemometer  also  failed  to  meet  sensitivity  requirements.  Analysis  of 
a  propeller  anemometer  using  blade  element  theory  made  clear  that  the  time  response  of  the 
instrument  will  be  proportional  to  the  product  of  density  and  velocity.  The  HAA  environment 
suffers  in  the  low  values  of  both  these  parameters.  A  propeller  with  very  acceptable  response  at 
low  altitudes  loses  approximately  an  order  of  magnitude  in  response  at  HAA  altitude  density  and 
speeds.  Winged  vehicles  that  maintain  nearly  constant  indicated  airspeed  with  altitude  would  be 
a  better  platform  for  pitot-static  and  propeller  anemometers. 

The  one  anemometer  that  appears  to  perform  better  at  higher  altitudes  is  the  ionic 
anemometer.  An  instrument  based  on  earlier  designs  of  French  investigators  was  built  and  tested 
in  our  laboratory.  The  5  to  6kV  potentials  required  to  achieve  a  good  ionic  flow  at  near  sea  level 
altitudes  was  reduced  to  a  bit  over  IkV  at  HAA  altitudes.  In  sea  level  calibration  testing,  the 
difference  divided  by  the  sum  of  the  ionic  currents  to  the  leeward  and  windward  collectors 
( Ai/Si )  proved  to  be  quite  linear  with  velocity  up  to  the  limits  of  the  wind  tunnel,  about  6  m/s. 

The  decision  to  build  an  inhouse  low  density  wind  tunnel  to  calibrate  the  ionic  anemometer  at 
HAA  operational  altitudes  proved  more  of  a  challenge  then  originally  expected.  Problems  were 
compounded  by  high  priority  activities  in  the  AFRL  shop  and  laboratories.  Consequently 
fabrication  of  the  device  was  delayed,  and  the  time  available  for  the  test  was  minimal.  In  spite  of 
this,  the  tunnel  worked  nearly  as  well  as  expected.  The  main  problem  was  leakage,  which 
prevented  operation  at  the  simulated  design  altitude  of  60  000  to  65  000  ft.  A  density 
representative  of  53  000  ft  was  sustained  long  enough  for  the  test.  Time  series  of  velocity  data 
were  collected,  and  a  representative  series  was  analyzed.  The  information  gained  by  statistical, 
PSD,  and  structure  function  analysis  of  the  data  showed  that  the  velocity  turbulence  data  was 
consistent  for  the  conditions  of  the  test,  and  appeared  to  have  the  levels  of  time  sensitivity 
attributed  to  ionic  anemometers,  which  is  the  detection  of  eddy  sizes  down  to  the  size  of  the  gap 
between  the  emitter  and  collector  electrodes. 
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The  ionic  anemometer  instrument  appears  quite  capable  of  making  velocity  fluctuation 
measurements  at  HAA  altitudes.  Actual  performance  has  only  been  verified  up  to  53  000  ft,  and 
a  using  agency  might  demand  a  simulation  up  to  full  design  altitude.  If  higher  altitudes  are 
required,  a  different  tunnel  must  be  used,  since  this  appears  to  be  the  low  pressure  limit  for  this 
PVC  pipe  material.  Potential  users  might  want  to  revisit  the  method  used  by  Barat,  which  was  a 
variant  of  the  first  configuration  analyzed  in  the  alternatives,  if  a  suitable  vacuum  chamber  were 
available  as  a  pressure  source.  In  addition,  a  notch  filter  should  be  used  around  the  60  Hz 
laboratory  noise  source  so  that  higher  frequency  information  such  as  the  inner  scale  can  be 
investigated  in  the  simulation. 
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Appendix  A:  High  Altitude  Wind  Tunnel  Start-up  Instructions  and 
Instrumentation  List 


The  following  outlines  the  instrumentation  and  general  proeedure  required  for  the 
sueeessful  start-up  and  operation  of  the  loeally  fabrieated  high  altitude  wind  tunnel  at  AFRL 
Hanseom  AFB,  MA. 

It  is  highly  reeommended  that  all  instrumentation  be  validated  as  operational  and  aeeurate 
utilizing  ambient  atmospherie  pressure  prior  to  initiating  a  vacuum  source  (to  arrive  at  desired 
altitude). 

Table  Al:  Wind  tunnel  Monitoring  Instruments  (for  locations,  see  Figure  Al) 


Nomenclature 

Method 

Measurement 

Range 

API 

Pitot  Tube  /  DAQ 

Velocity  via  AP 

0-1.0  mH20 

AP2 

DAQ  /  Manual 

Tunnel  AP 

0-2.0  mH20 

PI 

DAQ  /  Manual 

Tunnel  Mean  P 

0-760  mmHgA 

P2 

DAQ  /  Manual 

Vacuum  Pump  P 

0-760  mmHgA 

Tl 

DAQ  /  Manual 

Tunnel  Mean  T 

50-150  Deg  F 

RPM 

DAQ 

Fan  Speed 

As  Required 

VOLT 

DAQ 

Fan  Voltage 

As  Required 

AMP 

DAQ 

Fan  Current 

As  Required 

Atmospheric  Checkout 

Initial  check  out  of  tunnel  operation  and  instrumentation  should  be  accomplished  at  ambient 
sea  level  conditions  if  possible. 

1 .  With  the  tunnel  fully  instrumented  and  data  acquisition  system  active,  static  values  of  all 
instrumentation  should  be  checked  for  operation  against  known  zero  flow  calibration 
conditions  (that  is  AP’s  =  0,  T1  =  Ambient,  PI  =  Ambient,  VOLT  =  0,  RPM=0  etc.). 

2.  Fan  shall  be  placed  in  operation  and  rotational  velocity  slowly  increased  to  2,800  RPM. 
At  this  point,  tunnel  values  should  be  at  or  near  the  values  provided  in  Table  A2. 

3.  If  measured  values  are  not  in  agreement  with  Table  A2,  system  should  be  investigated  to 
determine  the  cause. 

4.  Closely  monitor  tunnel  temperature  Tl.  If  T1  exceeds  150  °F,  fan  shall  be  shut  down 
until  system  cools  to  acceptable  levels. 
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Table  A2:  Initial  Conditions  for  Sea  Level  Checkout 


Nomenclature 

Raw  Value 

Calculated  Value 

Ultimate  Calculated  Result 

API 

0.57  inH20 

50  ft/s 

Velocity 

AP2 

1.31  inHiO 

1.31  inH20 

Tunnel  Eoss  Efficiency 

PI 

'-'760  mmHgA 

-'760  mmHgA 

Air  Density  (altitude) 

P2 

-'760  mmHgA 

-'760  mmHgA 

Air  In-leakage  Validation 

Tl 

~70  °P 

~70  °P 

Air  Density  (altitude) 

RPM 

2,833 

1,055  EtVmin 

Volumetric  Elow  (velocity) 

VOET 

TBD 

270  Watts 

Tunnel  Input  Power 

AMP 

TBD 

270  Watts 

Tunnel  Input  Power 

Note: 

•  Due  to  the  much  higher  Power  requirements  of  the  fan  motor  at  high  density  (low 
altitude)  conditions,  heat  rejection  to  the  tunnel  will  be  several  times  greater  than  during 
low  density  (high  altitude)  operation.  Therefore,  tunnel  mean  temperature  Tl,  should  be 
closely  monitored  and  controlled  (using  cooling  air  bleed  if  needed).  Should  this  value 
exceed  150  °F,  tunnel  should  be  shut  down  and  allowed  to  cool  sufficiently  prior  to 
further  use. 

•  Due  to  the  very  low  AP’s  involved,  VERY  accurate  measurements  will  be  required  for 
these  values.  Recommend  angled  water  manometer  for  a  manual  cross  check  of  tunnel 
AP’s  and  pitot  tube  measurements. 


Initial  Checkout  at  Altitude 


Once  system  has  been  validated  under  sea  level  conditions,  tunnel  can  be  slowly  brought  to 
high  altitude  conditions  as  follows: 

1 .  Ensure  fan  is  not  in  motion. 

2.  Ensure  valve  VI  is  closed  and  that  valve  V2  is  fully  open. 

3.  Start  vacuum  pump. 

4.  Slowly  open  tunnel  isolation  valve  VI  while  monitoring  tunnel  pressure  PI. 

5.  Slowly  begin  to  close  pump  bleed  valve  V2  while  monitoring  tunnel  pressure  PI . 

6.  Continue  to  slowly  close  V2  until  the  desired  pressure  is  obtained  (~74  mmHgA) 

7.  Adjust  V2  to  maintain  optimum  pressure  throughout  the  run  cycle. 

8.  If  tunnel  pressure  remains  above  74  mmHgA  with  V2  fully  closed,  investigate  system  for 
air  in-leakage  and  correct.  Note:  Pump  performance  curve  (Eigure  A5)  can  be  used  to 
determine  the  actual  magnitude  of  air  in-leakage. 

9.  Start  fan  and  adjust  rotational  velocity  to  3,500  RPM 

10.  Adjust  fan  RPM  to  obtain  desired  flow  velocity.  If  flow  velocity  is  not  attained  by  4,500 
RPM  discontinue  test  and  increase  blade  angle  by  10°. 
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11.  Compare  measured  values  to  those  provided  in  Table  A3.  Signifieant  deviations  should 
be  noted  and  investigated. 

12.  Closely  monitor  tunnel  temperature  Tl.  If  T1  inereases  beyond  150  ‘’F,  an  artifieial  air 
leak  should  be  introdueed  just  downstream  of  the  8”  test  seetion  (by  opening  and 
adjusting  valve  V3)  to  cool  the  main  flow. 

Table  A3:  Initial  Conditions  for  60  000  ft  simulation 


Nomenclature 

Raw  Value 

Calculated  Value 

Ultimate  Calculated  Result 

API 

0.06  inH20 

50  ft/s 

Velocity 

AP2 

0.21  inH20 

0.21  inH20 

Tunnel  Foss  Efficiency 

PI 

74  mmHgA 

74  mmHgA 

Air  Density  (altitude) 

P2 

<74  mmHgA 

<  74  mmHgA 

Air  In-leakage  Validation 

Tl 

70  °F 

70  °F 

Air  Density  (altitude) 

RPM 

3,500 

1055  FtVmin 

Volumetric  Flow  (velocity) 

VOFT 

TBD 

80  Watts 

Tunnel  Input  Power 

AMP 

TBD 

80  Watts 

Tunnel  Input  Power 

Note: 

•  As  tunnel  pressure  decreases  it  is  important  to  closely  monitor  fan  RPM  as  an 
inadvertent  overspeed  due  to  the  lowering  density  conditions  could  cause  damage  to  the 
blades,  tunnel  walls  or  drive  motor. 

•  Due  to  the  very  low  densities  involved,  VERY  accurate  pitot  tube  measurements  will  be 
required  for  a  velocity  calculation.  Recommend  angled  water  manometer  for  a  manual 
cross  check  of  tunnel  AP’s. 

•  The  above  set  pressure  of  74  mmHgA  assumes  a  mean  air  temperature  of  70  °F.  The 
required  pressure  will  increase  as  cell  temperature  increases.  See  Section  4.3  in  body  of 
technical  report  for  discussion  and  plot  of  required  pressure. 

•  Final  fan  blade  angle  and  RPM  will  require  iterative  testing  once  system  has  been  check¬ 
out  on-line  at  altitude.  If  unable  to  obtain  required  volumetric  flow(velocity)  increase 
RPM  (up  to  a  maximum  of  4,500).  If  still  unable  to  generate  enough  flow  increase  blade 
angle  in  increments  of  1 0°  until  flow  is  obtained 
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Figure  Al.  Schematic  of  wind  tunnel  showing  locations  of  instruments  described  in  Tables  Al, 
A2,  and  A3. 
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Note:  Use  Teflon  Tape  on  all  threaded  connections. 
Once  installed,  lay  a  bead  of  Dow  Coming  High 
Vacuum^  grease  on  joints  if  high  air  leakage  is  noted 


1.5"  Brass  or  Stainless 
full  port  Ball  Valve  (VI) 


1.5  "  NPT 
Connection 


12" 

PVC 

Pipe 


Metered  Air  Leak  to 
maintain  vessel 
pressure  (0  to  30  Ib/hr) 
-OR-  just  install  a 
cheap  globe  \«l\e  and 
meter/throttle  as 
needed 


V 


O 


Discharge  to  Atmosphere 
(Optional  Air  Leakage  Meter 
can  be  added  here  as  well 
(but  not  used  during 
evacuation  ofvassel) 


Figure  A2.  Vacuum  pump  installation  detail 
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Fan  Code:  01 2/1 50/5/10 


3.0 


Requirements 

Volume: 

Static  Pressure: 
Selection  Pressure: 
Installation  Category: 
Temperature: 
Altitude: 


Fan  Data  (at  STP) 

Type 

Diameter: 

Hub: 

Impeller  Blades: 
Pitch: 

Blade  Material: 
Speed: 

Absorbed  Power: 
Peak  Power: 

Total  Efficiency: 


Motor  Data  (at  STP) 
Motor  Type: 

Electrical  Supply: 
Motor  Frame/Power: 
Current  FLC/Start: 
Motor  Speed: 


1000  cfm 
1  "WG 

1.12  "WG  at  std  conditions 
TYPE  D 


70  deg  F 
Oft 

Static 

Pressure 

CWG) 


Axial  Impeller 
12  in 
150  mm 
5 

10  degrees 
Aluminium 
3500  RPM 
0.34  hp 
0.36  hp 
56% 


Power 

(hp) 


II!' 

1  1  1  1 

1  1  1  1 

1  1  1  1 

1  1  1  1 

1  1  1  1 

1  1  1  1 
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Note:  Levels  are  quoted  as  in-duct  values.  dBA  values  are  average  spherical  free-field  for  comparitive  use  only. 
O  2001  Elta  Fans 

Figure  A3.  ELTA  fan  data  for  model  012/150/5/10 


Sample  ealeulation  of  required  RPM,  to  produee  a  desired  pressure  differenee; 

Fan  Law  Equations  (Known:  20°  blade  angle,  3500  RPM,  1055  ACFM,  2.0  inH20) 

CFM2  =  CFMl  X  (  RPM2  /  RPMl) 

AP2=APi  x(RPM2/RPMi)^ 

Power2  =  Poweri  x  (  RPM2  /  RPMi)^ 

Sea  Level  Example:  20°  blade  angle,  A  RPM,  1055  ACFM,  1.31  inH20 
1.31  =2.0  x(RPM2/3500f 

RPM7  =  2,832 
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Figure  A4.  Vacuum  pump  performance  curves  for  RA  and  RB  models. 


59 


Inle^  Presure-  Inches  Hg 


10  100  7eo 

Inlet  Pressure  -  Ton" 

Figure  A5.  Vacuum  pump  performance  curves  for  RC  models. 
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APPENDIX  B:  IONIC  ANEMOMETER  DATA  COLLECTION  ERROR 
ASSESSMENT 

Bl.  Introduction 

The  identification  and  measurement  of  clear  air  turbulence  in  the  higher  levels  of  the  atmosphere  (~  60 
000  ft)  is  of  prime  importance  for  the  USAF.  One  instrument  that  shows  great  promise  to  measure 
turbulence  data  (time  dependent  velocity  information)  is  the  Ionic  Anemometer. 

The  ionic  anemometer  is  a  unique  instrument  that  uses  a  high  voltage  source  to  promote  ion  movement 
between  charged  rods  (emitter  and  collector).  The  movement  of  these  ions  in  air  is  sensitive  to 
background  wind  velocity,  air  density  and  humidity  levels.  Activation  voltage  (the  minimum  voltage 
required  to  initiate  ion  motion)  is  strongly  a  function  of  air  density.  Lowering  the  air  density  also  lowers 
activation  voltage.  A  complete  theory  of  the  anemometer  operation  will  not  be  provided  here,  but  is 
available  in  the  literature 

The  ionic  anemometer  has  existed  at  the  research  level  for  several  decades  but  has  never  been 
developed  into  a  commercial  product.  Therefore,  attempts  to  deploy  of  this  type  of  instrument  typically 
require  local  manufacture  (based  upon  design  criteria  outlined  in  the  available  literature).  Very  limited 
data  exists  for  this  type  of  instrument  at  high  altitude  (low  density)  conditions.  As  this  instrument  showed 
great  promise,  VSBYA  fabricated  an  ionic  anemometer  in  2006.  This  device  was  tested  at  sea  level 
flowing  air  conditions  by  VSBYA  (see  additional  reports  in  this  TR).  Instrument  response  under  these 
conditions  was  excellent.  The  instrument  was  next  placed  in  a  vacuum  Bell  jar  to  simulate  response  at 
high  altitudes.  The  Bell  jar  test  was  able  to  establish  an  activation  voltage  vs.  pressure  curve  for  the 
instrument.  The  test  results  confirmed  that  activation  voltage  decreased  significantly  with  decreasing  air 
density.  While  in  the  Bell  jar  under  vacuum,  some  tests  of  the  instrument’s  ability  to  measure  air  velocity 
were  conducted  using  a  small  fan.  Though  the  arrangement  did  not  allow  for  a  quantitative  evaluation, 
the  results  did  prove  that  instrument  response  to  changes  in  air  velocity  was  quite  good. 

VSBYA  was  provided  an  opportunity  to  conduct  high  altitude  turbulence  data  collection  on  the  High 
Altitude  Airship  (HAA).  To  this  end,  a  better  understanding  of  the  instruments  response  to  velocity  data 
at  high  altitude  was  needed  to  better  quantify  the  instrument’s  suitability. 

In  2006,  visiting  intern  Mark  Duell  constructed  a  closed  loop  wind  tunnel  consisting  of  runs  of  12” 
PVC  piping  and  a  section  of  8”  PVC  piping  in  the  actual  test  section.  The  air  movement  was  provided  by 
a  12”  fan  with  a  variable  speed  drive.  Tunnel  pressure  was  set  by  an  externally  controlled  30  SCFM 
vacuum  pump  that  could  vary  vacuum  level.  Complete  details  are  available  in  the  body  of  this  technical 
report .  Data  were  collected  by  VSBYA  Captain  Barbeau  at  a  density  altitude  corresponding  to  138 
millibars  @  1 5  °C  and  also  at  sea  level,  the  results  of  which  are  reported  in  the  main  body  of  this 
Technical  Report. 

This  document  is  provided  as  a  means  to  estimate  uncertainty  associated  data  collected  while  using  the 
VSBYA  wind  tunnel 

B2.  Uncertainty  Overview 

Ultimately,  when  collecting  and  post  processing  data  the  question  that  arises  is  “how 
accurate  is  the  result”?  A  typical  approach  to  answer  this  question  is  to  conduct  an  uncertainty 
analysis.  Several  standards  are  available  to  conduct  such  an  analysis.  One  such  top  level 
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standard  is  International  Organization  for  Standardization  (ISO)  “Guide  to  the  Expression  of 
Uncertainty  in  Measurement” .  This  standard  was  published  in  1993  and  establishes  general 
rules  for  evaluating  and  expressing  uneertainty  in  measurements.  One  other  standard  often  used 
in  industry  is  AMSE  PTC  19. 1- 1998  “Test  Uncertainty” .  Though  the  ASME  standard  elosely 
follows  the  ISO  standard,  it  is  often  preferred  due  to  its  more  simplistie  formulations  and 
descriptions. 

Detailed  uncertainty  theory  will  not  be  covered  in  this  document.  The  reader  is  directed 
to  the  above  standards  for  further  information. 

In  general,  uncertainty  is  an  assessment  of  the  errors  associated  with  a  given 
measurement  (or  product  of  measurements).  Errors  in  measurement  arise  from  several  sources 
such  as  calibration  errors,  instrumentation  placement  errors,  data  acquisition  hardware/software 
interface  calibration  or  explicit  errors,  data  collection  methods  and  improperly  collected  transient 
information. 

These  errors  are  typically  divided  into  two  general  categories: 

1)  Systematic  Error  {a)  (often  referred  to  as  Bias  error)  is  a  portion  of  the  total  error  that 
will  essentially  remain  constant  for  repeated  measurements  of  a  parameter.  One  example 
of  a  systematic  error  would  be  an  instrument  calibration  error. 

2)  Random  Error  {/3  )  (often  referred  to  as  Precision  error)  is  a  portion  of  the  total  error  that 
produces  slight  variation  for  repeated  measurements  of  a  parameter.  These  variations 
will  show  a  distribution  of  data  for  a  collected  sample.  Given  enough  samples,  statistical 
methods  can  be  used  to  determine  an  assumed  true  value  of  the  measurement.  The  most 
common  distribution  of  random  error  is  the  Gaussian  or  normal  distribution.  Successive 
readings  will  cluster  about  a  central  value  and  will  extend  over  a  limited  interval 
surrounding  the  central  value.  Eor  a  Gaussian  distribution,  the  average  value  of  random 
error  is  zero.  Quite  often  a  Gaussian  distribution  is  assumed  when  evaluating  random 
error;  however,  for  this  assumption  to  be  used,  a  minimum  number  of  samples  (usually 
30  data  samples  or  better)  and/or  other  criteria  may  be  required  depending  on  type  and 
methods  of  data  collection. 

Total  error  (uncertainty)  combines  all  effects  of  Systematic  and  Random  error.  One  common 
way  to  calculate  uncertainty  is  through  the  method  of  propagation  of  errors.  This  method 
involves  establishing  sensitivity  factors.  Sensitivity  is  defined  as  the  response  of  the  final  result 
to  a  slight  variation  in  one  of  the  measured  parameters.  Sensitivity  factors  can  be  calculated 
numerically,  analytically,  or  in  some  cases  measured  directly  from  experimental  data.  Eor  an 
example;  if  our  end  result  G,  relied  on  (or  was  a  function  of)  measurements  of  parameters  x,  y, 
and  z;  then  the  sensitivities  to  parameters  x,  y,  and  z  would  be  as  follows: 

Assume  G  (x,y,z)  =  (5x)(3y)/(z) 


Sensitivity  of  x  on  G  =  AG/ Ax  (at  small  Ax)  or 


dx 


15t 
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Sensitivity  of  y  on  G  =  AG/ Ay  (at  small  Ay)  or  —  =  — ^ 

dy  z 


Sensitivity  of  z  on  G  =  AG/Az  (at  small  Az)  or 

dz  z 

As  one  might  expeet,  measured  parameters  that  produce  the  largest  sensitivities  produce  the 
largest  uncertainties  in  the  end  result. 

Combining  the  effects  of  uncertainty  is  not  just  a  simple  matter  of  addition  (see  previous 
referenced  standards).  For  a  sufficient  number  of  samples  (>  '-SO)  the  following  procedure  is 
often  employed  to  indicate  95%  confidence: 

Uncertainty  of  an  individual  measured  parameter  “x”: 

Where  a  (systematic  error)  and  P  (random  error)  are  known  (or  estimated)  based  upon 
calibration  information  and  statistical  sampling. 

Uncertainty  of  a  calculated  result  (of  more  than  one  measured  parameter): 


Utotal 


+  ...  or  X 


di 


The  reporting  of  a  calculated  value  (end  result)  will  include  the  calculated  value,  the  uncertainty 
and  the  confidence  level  used  to  obtain  the  uncertainty  level.  For  example  the  measured  value  of 
heat  load  to  air  (Btu/hr)  may  have  uncertainty  expressed  as  a  direct  value  or  as  a  percentage 
value: 


Heat  Load  “Q”  =  Qmeasured  (Btu/hr)  +/-  110  Btu/hr  (95%) 

or 


Heat  Load  “Q”  =  Qmeasured  (Btu/hr)  +/-  4.5%  (95%) 


B3.  Governing  Equations  and  Measured  Parameters  for  the  Ionic 
Anemometer 

The  ionic  anemometer  is  being  employed  to  measure  velocity.  The  value  of  velocity  is  to 
be  inferred  from  voltage  and  current  measurements  to  and  from  the  instrument.  Corrections 
must  also  be  applied  due  to  changes  in  density,  applied  voltage  and  humidity.  Ultimately  mean 
velocity  and  mean  ionic  anemometer  current  levels  are  evaluated  to  obtain  the  proportionality 
constant  K,  which  will  complete  the  equation: 
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V  =  K(l2-Il)/(l2+Il)  +  B 


where  v  is  veloeity.  Note  that  K  and  B  will  have  units  of  veloeity.  B  will  be  an  experimentally 
determined  zero  offset,  /y  and  I2  are  upwind  and  downwind  eolleetor  eurrents  respeetively. 

The  VSBYA  lonie  Anemometer  used  for  the  wind  tunnel  evaluation  eonsisted  of  a 
speeialized  eurrent  sensing  paekage  developed  by  VSB  eleetrieal  engineer  Kris  Robinson.  This 
sensing  paekage  reported  (I2-I1)  =  A  and  (I2+I1)  =  E,  direetly,  in  lieu  of  Ii  and  I2  individually. 
This  was  done  in  an  effort  to  inerease  aeeuraey  of  the  end  result  and  aid  in  flexibility  of  data 
eolleetion.  This  deviee  also  allowed  for  pre-zeroing  of  output  data  (in  a  known  zero  veloeity 
flow),  henee  removing  our  offset  eonstant  B  from  the  equation. 

Our  new  equation  beeomes: 


In  general  our  overall  assessment  will  be  the  result  of  tunnel  inputs  from: 

•  Instrument  Input  Voltage  (V) 

•  Instrument  Colleetor  Current  Differenee  (A) 

•  Instrument  Colleetor  Current  Sum  (E) 

•  Tunnel  Absolute  Pressure  (P) 

•  Tunnel  Fan  Speed  (co),  used  to  direetly  indieate  tunnel  wind  Velocity 

•  Tunnel  Temperature  (T) 

Additional  errors  in  final  result  from: 


•  Plaeement  of  pressure  taps 

•  Unknown  tunnel  turbulenee  levels 

•  Assumptions  of  aeouraey  with  fan  law  equations  in  a  low  density  environment 

•  Single  point  data  sampling 

•  Un-ealibrated  pressure  instruments 

•  Unknown  tunnel  humidity  levels 

•  Tunnel  Air  Leaks 

•  Impreeise  tunnel  pressure  eontrol 

•  Non-exaet  instrument  orientation  in  tunnel  (ionie  anemometer) 

•  Unknown  tunnel  veloeity  profile 

We  also  know  that  A  and  E  are  affeeted  by  ion  mobility.  Ion  mobility  is  strongly 
infiueneed  by  density,  humidity  and  applied  voltage.  Therefore,  all  variables  must  be  taken  into 
aeeount  in  the  uneertainty  analysis. 


Governing  equation(s): 
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V 


(A)/(Z) 

Where; 

K  =  Proportionality  Constant 

V  =  Mean  tunnel  wind  speed  — >  F(co  only) 

A  =  Average  Colleetor  Current  Difference  — >  F(P,  T,  Voltage(V)  and 
Humidity(H)) 

E  =  Average  Collector  Current  Sum  — >  F(P,  T,  V,  and  H) 

The  above  implies  that  K  =  K(V,  P,  T,  A,  E,  co,  H) 

Overall  uncertainty  will  take  the  form: 


Utotal 
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B4.  Estimating  the  individual  Ui’s 


Values  for  Ui’s  were  computed  utilizing  OEM  values  for  instrument  accuracies  and 
general  estimates  from  the  users  (for  a  and  P  values).  The  following  have  been  used  as  a  basis 
for  this  evaluation: 


Table  B1 :  Uncertainty  estimates  for  velocity  measurements 

a  Uncertainty 

Value  Units  Source 

P  Uncertainty 
Value  Units 

Source 

U„a=  1.08 

%  Reading  Vendor  Data 

UvP  = 

500 

Volts 

Estimate 

Upa=  5.0 

millibar  Estimate 

UpP  = 

10.0 

millibar 

Estimate 

UTa=  0.5 

°C  Estimate 

UtP  = 

2 

°C 

Estimate 

U^a=  1.05 

%  Reading  Vendor  Data 

UaP  = 

0.03 

mA 

Estimate 

Uj;a=  1.05 

%  Reading  Vendor  Data 

U2P  = 

0.05 

loA 

Estimate 

U„a=  1.05 

%  Reading  Vendor  Data 

U„P  = 

1 

Hz 

Estimate 

Uh  a  =  50 

%  Estimate 

UhP  = 

50 

% 

Estimate 

The  ultimate  values  for  Ui  =  (ai^  +  Pi^)°  ^  will  typically  be  unique  for  each  data  point 
when  accuracy  values  are  provided  in  %  of  reading,  hence  a  Ui  table  is  not  provided  here. 


B5.  Establishing  Sensitivity  Factors  {dKjdxp 

Evaluating  our  governing  equation  for  all  the  variables  listed  above  K(co,A,E,V,  P,T,H)  is 
accomplished  using  a  combination  of  analytical  and  empirical  methods.  The  results  are  provided 
here,  with  details  provided  below. 

Tunnel  initial  test  calibration  runs  readily  arrived  at  an  expression  for  Velocity  in  terms  of 
fan  motor  frequency  co.  Though  this  probably  could  have  been  arrived  at  in  an  analytical 
fashion,  the  direct  calibration  of  the  tunnel  against  a  precision  turbine  wind  meter  was  a  simple 
expedient  and  provided  the  following  linear  relationship.  Note,  a  linear  relationship  was  also 
expected  due  to  incompressible  fan  theory  and  should  also  remain  valid  for  lower  density 
conditions). 

Velocity  (m/s)  =  0.0206439(m/(s-Hz))*co(Hz) 

Our  governing  equation  can  now  be  written  as: 

0.020643907 

(A)/(s) 
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For  the  explicit  variables  co,  A  and  E,  the  needed  partial  derivatives  can  readily  be 
evaluated  as: 


—  =  0.0206439  %  (m/s(Hz)) 

dm  /  ^ 


—  =  -0.0206439^%  (m/(s-pA)) 

DA  /A 


dK 


-0.0206439  ^ 


(m/(s-|uA)) 


Analytical  expressions  for  dependence  of  K  on  V,  P,  T  and  H  are  not  readily  available. 
However,  research  from  previous  authors  ’  ’  ,  and  our  test  data  can  provide  meaningful  insight 
when  estimating  sensitivity  factors  for  these  variables.  Small  perturbations  in  actual  test  data  for 
the  variable  in  question  (with  all  other  variables  held  constant)  can  also  provide  a  good  estimate 

r-  1  ,dK  AK  ^ 

ot  the  needed  derivatives  ( —  = - ). 

AX; 

A  review  of  J.  Barat’s  paper  A  High  Resolution  Anemometer  for  Boundary-Layer 
Measurements  reveals  sensitivity  assessments  based  upon  measured  response  to  parameter 
changes.  Section  4  of  his  paper  suggests  the  following: 

K  sensitivity  to  Voltage  Changes: 


K  =  Ky\  +  l.iE-^V-V,)] 


Where: 

Ko  =  Constant  corresponding  to  V=Vo 

Vo  =  5700Volts  (for  sea  level  conditions  and  his  geometry) 

This  formulation  was  compared  to  our  sea-level  test  data.  Following  Barat’s 
methodology,  an  adjustment  was  made  to  the  voltage  constant  yielding: 

A  =  Ao[l  +  3.8£-4(F-Fo)] 
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As  we  were  unable  to  collect  good  voltage  perturbation  data  at  altitude,  we  will  assume 
that  this  sea-level  sensitivity  holds  at  altitude.  This  sensitivity  is  provided  as; 

=  3.8x1 (m/(sV)) 
oV 

(Note:  VSB  Ko  value  is  30.70  m/s  at  sea-level) 

K  sensitivity  to  changes  in  Pressure  Changes: 

Note;  Barat  calls  this  sensitivity  to  density,  though  he  only  varies  pressure. 


K  =  K^ 


1  +  0.07 


Po-P 


)) 


Where; 


K'^  =  Constant  corresponding  to  P=Po 

Po  =  1013mb  (for  sea  level  conditions  and  his  geometry) 


This  would  suggest; 


dK  -K'O.OlPo  „  ,  , 

—  = - - — ; - m/(s  mbar) 

dP  P^ 


This  formulation  is  predicated  on  V  =  Vo  =  constant.  Unfortunately,  VSB  data  collection 
did  not  provide  a  ready  means  to  verify  this  formulation  (I.e.  data  not  collected  for  pressure 
perturbations  with  voltage  held  constant).  VSB  data  did,  however,  provide  pressure 
perturbations  with  voltage  set  at  activation  for  each  pressure.  Though  not  ideal  (due  to  the 
varying  voltage)  this  sensitivity  factor  has  been  estimated  using; 


dK 

dP 


AK 

AP 


=  -7.385x10^' 


(m/(s  mbar)) 


K  sensitivity  to  temperature: 

We  know  from  literature  that  ion  mobility  is  strongly  a  function  of  the  density  of  the 
flowing  medium.  The  literature  does  not  address  any  difference  in  density  caused  by 
temperature  changes  or  pressure  changes.  As  we  cannot  measure  density  directly  (it  becomes  a 
calculated  value  from  a  measured  P&T),  we  can  only  assume  that  the  sensitivity  found  for 
pressure  will  be  related  to  the  sensitivity  for  temperature  owing  to  an  appropriate  equation  of 
state.  As  we  are  dealing  with  air  at  ambient  temperatures  (~  15  °C),  the  Ideal  Gas  Law  has  been 
imposed  for  this  purpose  (P  =  pRT,  where  R  is  the  MW  adjusted  gas  constant  for  air).  Utilizing; 
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The  Ideal  Gas  Law  provides: 


dK  ^  dK  dP  dp 
dT  dP  dp  dT 


dK 

dT 


dp 

dT 


SP  r^rr. 

—  =  RT 
dp 

P 


Rpl 

-P 


—  {RT)  2 
dP  ^  ’  RT^ 


dK  -P 
dP  T 


For  the  conditions  of  our  test,  the  value  of  dKjdP  is: 


dK  P 

—  ^7.835x10^'-  (m/(sK)) 

dT  T 


K  sensitivity  to  humidity  changes: 

It  is  well  known  that  ion  mobility  is  sensitive  to  humidity  levels  (mobility  decreasing 
with  increasing  humidity  levels).  Barat  has  provided  the  following  based  upon  his  collected  test 
data  at  sea  level. 


=  (1-0.03577) 

Where: 


7f"  =  Constant  corresponding  to  77=  0.0 
77=  Relative  Humidity  from  0.2  to  0.7 


This  would  imply: 

dK  _ 
dH  ~ 


-.^;(0.035)  (m/s) 


for  50%  humidity.  For  our  use,  we  will  use  7f"  =  92.71  m/s  (our  measured  K  at  the  138  millibar 
test  conditions). 


Note:  The  above  is  formulated  on  Relative  Humidity  at  sea-level.  As  relative  humidity 
changes  with  pressure  and  temperature  (for  a  given  amount  of  moisture),  it  may  make  more 
sense  to  formulate  in  terms  of  Specific  Humidity. 
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B6.  Overall  Uncertainty  of  a  Test  Point 


The  above  uncertainties  and  sensitivities  can  now  provide  an  estimate  of  the  overall 
uncertainty  of  our  proportionality  constant  K  for  each  data  point.  This  value  can  then  be 
translated  to  an  uncertainty  in  a  calculated  velocity  by  using  the  relations: 


Velocity  =  —  (if  ±  )  or 


Velocity 


B7.  Results  at  the  138  mbar  test  condition 


The  above  formulations  have  been  utilized  with  our  low  pressure  test  data  at  138 
millibars.  The  plot  below  provides  the  uncertainty  estimates  as  error  bars  on  the  plot. 
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Figure  Bl.  Conversion  curve  of  mean  windspeed  vs.  the  mean  ratio  of  the  difference  of  current 
to  the  sum  of  the  currents  with  estimated  error  bars  on  the  data.  Also  shown  is  the  anticipated 
behavior  of  the  curve  based  on  published  performance  of  the  ionic  anemometer. 


The  plot  provides  a  dotted  line  corresponding  to  the  K  value  used  for  the  sensitivity 
factors  for  this  data  set.  From  previous  work,  it  is  fully  anticipated  that  the  instrument’s  velocity 
response  should  be  linear  with  respect  to  ratio  (down  to  and  including  the  zero  velocity  point). 
The  indicated  deviation  below  a  ratio  of  0.0375  will  require  further  review  to  give  us  a  better 
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understand  the  errors  assoeiated  with  these  partieular  data  points  (as  they  fall  well  outside  our 
estimated  uneertainty  bands).  Data  and  eorresponding  uneertainties  provided  for  ratio  values 
greater  than  0.0375  should  remain  valid  however. 

The  above  uneertainty  review  addresses  steady  state  mean  values  for  eolleeted  data  only, 
and  does  not  address  uneertainties  involved  when  using  the  instrument  for  turbulenee  data 
eolleetion.  It  is  assumed  that  the  uneertainty  assoeiated  with  the  mean  values  ean  also  be  applied 
to  the  turbulenee;  however  this  was  not  verified. 
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NOMENCLATURE 


a  = 

b 

c  = 

Cl 

Cu' 

D 

Duu(r)  = 
E 
Gr 
g 

i  = 

J 

K 

L 

M 

Nu 

Pr 

P,  p  = 

R 

r  = 

s  = 

T 

U,V  = 

V  = 

Z 

a  = 

A 

e 

+ 

p 

V  = 

p 

Z 

0 

Q 

Subscripts: 

b 

i  = 

M 

o  = 

P 

P 

S 

1,2 

00  = 


Lift  curve  slope  or  a  constant 
Number  of  blades 
Chord  (m) 

Coefficient  of  Lift  (n.d.) 

Structure  Constant  for  veloeity  (m'^^^)/s^ 

Drag  (N)  or  Diameter  (m) 

Strueture  Funetion  for  veloeity  over  a  distance  r  (m  /s  ) 

Eleetric  field  (V/m) 

Grashof  Number  (n.d.)  (defined  in  Eq.  27) 

Aceeleration  of  gravity  (m/s  ) 

Ion  current  (A) 

Moment  of  inertia  (kg  m  ) 

Thermal  conduetivity  (J/m  K) 

Eift  (N) 

Moment  on  the  propeller  shaft  (Nm) 

Nusselt  Number  (n.d.) 

Prandtl  Number  (n.d.) 

Pressure  (Pa  or  millibars  =  hPa) 

Radius  of  propeller  blades  (m) 

Radius  or  distance  (m) 

Distanee  along  a  path  (m) 

Torque  from  propeller  (Nm) 

Velocity  (m/s) 

Velocity  (mean  or  perturbation)  (m/s) 

Height  (m) 

Angle  of  attack  (deg)  or  thermal  diffusivity  (m  /s) 

Difference  of  the  two  ionic  collector  currents  (A) 

Pitch  angle  of  propeller  blade  (deg) 

Mobility  of  ions  (m^/sV) 

Kinematie  viseosity  (m^/s) 

Density  of  the  air  (kg/m^) 

Sum  of  the  two  ionic  collector  currents  (A) 

Twist  angle  (deg) 

Shaft  angular  veloeity  (rad/s) 

Velocity  entering  blade  element 
Index  for  two  ionie  veloeity  terms 
Motor 

Reduced  mobility 
Propeller 

Veloeity  entering  front  of  propeller  anemometer 
Shaft 

Colleetor  elements  of  ionic  anemometer,  or  loeations  along  a  streamline 
Baekground  flow  (veloeity) 


73 


